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FOREWORD 


The thermal analysis of Skylab fuel cell and 
radiator system capability, in support of Contract 
SA500, NAS9-150, was conducted between October 1969 
and June 1970. Results of this work are contained in 
this report, which was prepared by J. L. Schaefer 
and G.A. Vanderpol of CSM Operations Analysis, 
Mission Requirements and Evaluation group. The 
original models of the electrical power system (EPS) 
radiator and Block I fuel cell, and the CINDA com- 
puter program were furnished by H. Cazemier and 
W. Simon of NASA-MSC. Early checkout runs were 
limited to leased Univac 1108 service and IBM 7094 
emulation, which were slow and expensive. With the 
assistance of P, Jepsen of Aero and Thermal Projects, 
Scientific Programming, the IBM 360/SINDA program 
was obtained and converted to Space Division's IBM 
360 formats. E. R. Arnold provided considerable 
assistance in developing the Block II fuel cell model 
for the SINDA program. 
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1. INTRODUCTION 


1. 1 PURPOSE OF ANALYSIS 

The range of output power capability for the Apollo fuel cell system is 
established by the requirement to operate within the command module (CM) 
bus voltage range of 26- 5 to 31 volts and by the thermal constraints imposed 
by the fuel cell stack, condenser, and radiator temperature limits. This 
report considers only the Uiermal constraints in establishing the fuel cell 
power capability for Skylab missions. Under high transient -power levels, 

CM voltage requirements may be more limiting. 

Early in 1969, NASA-MSC personnel furnished for the Skylab study the 
computer program, the Chrysler Improved Numerical Differencing Analyzer 
for Third Generation computers (CINDA-3G) (Reference 1), and heat transfer 
models of the Block II electrical power system (EPS) and the Allis -Chalmers 
fuel cell system. Subsequently at NR several major changes were made to 
these models and the CINDA program. The preliminary Allis-Chalmers 
model became obsolete after the decision was made to use the Pratt and 
Whitney (PW) Block II fuel cell for the Skylab mission. Therefore, a com- 
pletely new PW fuel cell model had to be developed for CINDA. Several 
modifications were made to the EPS radiator model to provide closer corre- 
lation with qualification test data (Reference 2) and to include heat transfer 
between the radiator panels and the service module (SM) structure. Finally, 
the original CINDA- 3G program, which was developed for use on the 
Univac 1108, was replaced with an advanced program version, Systems 
Improved Numerical Differencing Analyzer for Third Generator Computers 
(SINDA-3G), which had been converted (Reference 3) for operation on the 
IBM 360/75 computer, 

1.2 VERIFICATION ANALYSIS REQUIREMENTS 

Verification analysis requirements consist of a computer simulation. 

The simulation determines the effect of the thermal constraints on power 
capability with one fuel cell and with two fuel cells and with full and five- 
eighths EPS radiator area operation under the two extremes for Skylab fixed 
orbit environmental conditions. The SINDA-3G computer program is required 
for the eight cases of computer simulation. Results for each of the eight 
cases must be within the acceptable fuel cell temperature limits stated for 
the flight measurements. 
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2. CONCLUSIONS 


Analysis results indicate the EPS radiator subsystem has sufficient 
capability to reject the waste heat associated with the fuel cell power levels 
required for Skylab missions while maintaining the fuel cells within nominal 
temperature levels. The computer simulation considered sun-vector orbit 
plane angles (p) of 0,0 and 73. 5 degrees for one and two fuel cells operating 
with full and five-eighths radiator areas. The maximum and minimum 
power levels evaluated in this study are listed in Table 1. Results are dis- 
cussed in Section 4. 4 of this report. 


Table 1* Sumnaary of Results for Fuel Cell/Radiator Power 
Capability for an Earth Orbit of 235 Nautical Miles 


Fuel Cells 
(number operating) 

Radiator Area 
(operating panels) 

Minimum 
Total Current 
(amperes) 

Maximum 
Toteil Current 
(amperes) 

. 2 

8 

50 , 

100 

2 

5 

40 

90 

1 

8 

30 

60 

1 

5 

25 

50 
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3. RECOMMENDATIONS 


Based upon results of this anetlysis, qualification and development test 
results, and Block II flight data, Block II fuel cell and EPS radiator sub- 
system capabilities ‘are sufficient to meet Skylab mission requirements. No 
further testing or analysis is recommended. 


Preceding page blank 
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4. ANALYSIS' 


4. I ANALYTICAL MODELS 

General Fuel Cell Model Considerations 


The Pratt and Whitney fuel cell combines oxygen and hydrogen to pro- 
duce electricity, heat, and water. The overall chemical reaction can be 
written as 


H 2 + 1/2 O 2 — » H 2 O + Energy 

The quantity of reactants required to produce a given amount of energy can 
be determined from Faraday's Laws of Electrolysis. As stated in Refer- 
ence 4, these are: 

1. The mass of a substance liberated in an electrolysis cell is pro- 
portional to the quantity of electricity passing through the cell. 

2. When the same quantity of electricity is passed through different 
cells, the masses of the substances liberated are proportional to 
their electrochemical equivalents. 

The cprnbined laws can be used to express the amount of water pro- 
duced per unit time in terms of current flow: 


Wh^O = I-M/(F=^n) (1) 

where 

I = current (amperes) 

M = molecular weight 

F = Faraday's constant— 96500 amp- sec /gm equivalent 
n = gm equivalent/gm mole 
W = weight flow per unit of time 

Preceding page blank_ 
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The following values of specific fuel consTimption have been determined 
from the above equation: 


0.0230 lb H2O per amp-hr 
0.0204 lb ©2 per amp-hr 
0. 00257 lb H2 per amp-hr 

Figure 1 depicts the basic schematic of the fuel cell thermodynamic 
components. Two coolant loops are used to remove the excess heat and 
water. The stack contains the 31 individual fuel cell elements, each of 
which consists of two electrodes and the KOH electrolyte, where the chemical 
reaction of hydrogen and oxygen occurs. Electrically the cells are connected 
in series. Each is capable of approximately 1 volt, depending on electrolyte 
condition and load. The temperature of the stack is maintained in the range 
of 390 to 460 F by means of the primary regenerator and bypass valve. The 
condenser and water separator are responsible for the water removal and the 
resulting water concentration in the electrolyte. The temperature range at 
which water is nominally condensed is 155 to 165 F. 

This temperature range is controlled by the secondary regenerator 
bypass valve and the secondary regenerator. After some of the waste heat 
is used to heat the incoming reactants in the oxygen and hydrogen preheaters, 
the EPS radiators reject the excess heat to space. The analytical modeling 
for each major component is discussed in more detail in the next sections. 

Voltage Output 


The fuel cell terminal voltage is a function of load current, stack 
temperature, and electrolyte (KOH) water concentration. Reference 5 
provides nominal fuel cell performance curves that depict voltage as a func- 
tion of current and parametric values of stack (surface) temperature. These 
curves were expressed in the general polynomial form 

Vt = Aj + A2 I + A3l^ + A4l^ + A5TS + A6(TS) (I) 

+ A7(TS)(I^) + AgTS^ + A9 (I)(Ts 2) f AiqTS^ (2) 


where 


Vj. is voltage based on current and stack temperature (volts) 

Aj Ajq are coefficients 

I = load current (amperes) 

TS = stack temperature (°F) 
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A standard least squares surface fit bivariate polynomial routine was 
used in calculating the coefficients of this equation. The Input data for the 
bivariate routine, as shown in Table Z, are the Reference 5 data corrected to 
a common KOH-water concentration of 0.75. The coefficients established 
from these data are given in Table 3. The resulting value, V^, is corrected 
for the KOH- water concentration effect in the SIN DA program by the following 
equation: 


V = -0.241 (0. 75 - PCKOH) 


(3) 


where 


V = fuel cell voltage based on current, stack temperature, and 
KOH -water concentration (volts) 

Vj. = voltage based on current and stack temperature (volts) 

PCKOH = the weight ratio of KOH to KOH + water 

The value of V is used In the program as the fuel cell terminal voltage. 


Table 2. Input Data for Bivariate Routine 


Current 

(amperes) 

Stack Temperature (“F) 

380 

400 

420 

440 

460 

15.0 

30.0 

30.8 

31.4 

31.8 

31.9 

20.0 

29. 1 

30.0 

30. 6 

31.0 

31.3 

25.0 

28.3 

29. 2 

29.9 

30.4 

30.7 

30.0 

27. 5 

28.5 

29.4 

29.8 

30. 1 

35. 0 

26; 7 

27. 8 

28. 8 

29.2 

29. 6 

40. 0 

26.0 

27.2 

28.2 

! 

28. 7 

29.0 

45.0 

25. 3 

26.5 

27.6 

28.2 

28. 5 

50. 0 

24.6 

25. 9 

27.0 

27. 7 

28. 1 
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Table 3. Output Coefficients From Bivariate Routine 


Aj 5.583364 

A^ 4.686233 x 10“3 

A 2 -1.308438 

A 7 -3.573639 x lO-"^ 

A 3 2, 054927 X lO"^ 

Ag -1.525995 x 10-5 

A 4 -1.333418 x 10-5 

A 9 -4. 889387 x 10~6 

A 5 9.694690 X 10 '^ 

A^q -1.292755 X lO""^ 


Primary Loop Thermodynamics 


Excess heat and steam produced in the chemical reaction in the stack 
are removed by the recirculating stream of hydrogen and water. The energy 
balance computation for the stack is based upon the total energy of reaction 
as input to the stack and the electrical power, stack temperature change, 
recirculating fluid temperature change, and fuel cell heat loss as output 
energy. The reaction energy is determined by 

QGENAT = 51600. SFCH2 (4) 


where 


QGENAT = reaction energy (Btu/hr) 

SFCH2 = hydrogen consumption rate (Ib/hr) 

51600 Btu/lb is the lower heating value of hydrogen 
The heat balance on the stack is given by 

QSTORD = QGENAT - QH2 - Q02 - QELECT - QRS - QSM (5) 


where 


QSTORD = heat gained in the stack (Btu/hr) 

QGENAT = reaction energy (Btu/hr) 

QH2 = heat gained by the consumed hydrogen (Btu/hr) 
Q02 = heat gained by the consumed oxygen (Btu/hr) 
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QELECT = electrical energy output (Btu/hr) 

QRS = heat gained by the recirculating hydrogen-water stream 
(Btu/hr) 

QSM = heat lost by the stack to the SM structure (Btu/hr) 

A similar equation can be expressed for the primary loop heat balance: 
QCOND = QGENAT - QELECT - QSM - QH20 (6) 

where 

QCOND s heat transferred across the condenser (Btu/hr) 

QH20 = energy in the condensed water above the 70 F datum 
(Btu/hr) 

QGENAT, QELECT, and QSM are as previously defined. 

After the two equations are balanced, the temperature can be deter- 
mined in the primary loop by the following three equations: 

TSE = TSTACK + (QSTORD=J‘DTIMEU/30. ) (7) 

where 

TSE = the new stack temperature { “F) 

TSTACK = the previous stack temperature ( “F) 

QSTORD = the stack heat gain (Btu/hr) 

DTIMEU = the time step (hr) 

30. = the stack mass -specific heat product (Btu/“F) 

TSI = TSE - QSR/(WDTCPI + WDTCP2) (8) 

where 

TSI = recirculating stream temperature at stack inlet ("F) 

TSE = recirculating stream temperature at stack outlet (®F) 

(as well as the stack temperature) 
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QSR = heat gained by the recirculating hydrogen-water stream 
(Btu/hr) 

WDTCPl = the mass -specific heat product of the water in the 
recirculating stream 

WDTCP2 = the mass -specific heat product of the hydrogen in the 
recirculating stream 

TCIP = TCEP - (QCOND - CONST2)/(WDTCPl h WDTCP2) (9) 

where 

■ 

TCIP = recirculating stream temperature at condenser primary 
inlet (*F) 

TCEP = recirculating stream temperature at condenser primary 
outlet { *F) 

C0NST2 = heat of vaporization of the condensed water 

QCOND, WDTCPl, WDTCP2 are as previously defined. 

The TCEP is determined in the condenser subroutine, which is 
described in the next section. 

. The mass balance in the primary loop is determined by assuming a 
constant volume delivery of 3 cfm at 60 psia for the primary pump for cal- 
culating the specific volume of the hydrogen and water at the stack inlet and 
outlet and at the condenser inlet and outlet. The partial pressure of water 
is first determined at the condenser exit, based upon the condenser exit 
temperature and a saturated steam condition. At the stack, the recirculating 
steam enters at this same condition and leaves at the temperature and the 
partial pressure of water in the electrolyte. Figure 2 illustrates the relation- 
ship of the electrolyte temperature, partial pressure of water, and KOH 
concentration. The program computes the equilibrium water pressure at 
the condition of electrolyte temperature and concentration at the beginning 
of the time step. This pressure is assumed to equal the partial pressure of 
water in the hydrogen -water stream at the stack exit at the end of the time 
step. Thus, the water balance is determined froth the following equations; 

WH202P = I>!=0. 0230 DTIMEU (10) 
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TEMPERATURE (»F) 

Figure Z. Water Pressure of KOH Electrolyte 
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where 

WH202P = weight of water produced within the time step (lb) 

I = current (amps) 

DTIMEU = time step (hours) 

0. 0230 is the-conversion factor for pounds of water produced per 
ampere-hour 

WH2025 = WH2026 + WH202P RATE , (11) 

where 


■WH2025 = weight of water in the electrolyte at the end of the time 
step (lb) 

WH2026 = weight of water in the electrolyte at the start of the 
time step (lb) 

RATE = the percent of water produced during the time step that 

remains in the stack. RATE is calculated by an iterative 
balance between the water specific volume at the stack 
inlet and at the stack outlet and the water production rate. 

PCKOH = 22. 0/(22. 0 + WH2025) (12) 


where 


PCKOH = electrolyte concentration (ratio of KOH weight to total 
electrolyte weight) 

KOH weight is 22. 0 pounds 

WH2025 is as previously defined 

The mass balance of water at the condenser uses the specific volume 
of water determined at the stack outlet as the specific volume of water at the 
condenser inlet. The specific volume of water at the condenser outlet is 
determined by the condenser primary exit temperature, TCEP, and a 
saturated steam condition. 
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The amount of water condensed is found as follows; 


DMASS = (180. /SVOLll)*DTIMEU - (180. /SVOL21)^'=DTIMEU (13) 


where 


DMASS = water condensed (lb) 

SVOLll = specific volume of water at the condenser outlet 
(cu ft/ lb) 

SVOL12 = specific volume of water at the condenser inlet 
(cu ft/ lb) 

DTIMEU = time step (hours) 

The mass balance for the reactants is based on providing the necessary 
flow to support the theoretical consumption rates without consideration of 
minor variations in regulated pressure: 

SFCH2 = 8,292 ^ 10"^ >!= 31. ^ I (14) 

SFC02 = 7.94 »= SFCH2 (15) 


where 


SFCH2 = hydrogen specific fuel consumption (Ib/hr) 

SFC02 = oxygen specific fuel consumption (lb /hr) 

The energy balance involving the consumed reactants is discussed in 
the Reactant Preheater section. 

Condenser 


The condenser is a counterflow heat exchanger interfacing the primary 
and secondary loops of the fuel cells. The condenser transfers the heat of 
the hydrogen-water mixture of the primary loop to the water -ethylene glycol 
rrdxture of the secondary loop. The normal operating range of the condenser 
exit temperature on the primary side (TCEP), which controls the secondary 
regenerator bypass valve, is 155 to 165 F. 

Several schemes were considered for defining condenser performance. 
The method used, in view of simplicity of interface between the primary and 
secondary loops, was the balance of heat flow across the condenser. Data 
describing the condenser were obtained from Reference 6 and are provided in 
Figure 3 in the form of condenser performance curves. A multiple linear 
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Figure 3, Condenser Performance Curves 
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regression scheme was then used to curve-fit the data. The following 
equation was obtained: 

TCEP-TCES = (1. /WDWG)*(49.229 + 0. 68747 QCOND +9.0145E-7 

QCOND ^ ^ Z) (16) 


where 


TCEP = condenser primary side exhaust temperature ( "F) 

TCES = condenser secondary side exhaust temperature (®F) 

WDWG = water -glycol flow rate (Ib/hr) 

QCOND = heat flow across the condenser (Btu/hr) 

The program makes successive approximations to balance the heat 
flow across the condenser. A value of QCOND is calculated from the primary 
loop parameters. Then the condenser glycol and gas exit temperatures are 
varied in accordance with Equation (16) until the heat lost to the glycol equals 
the QCOND calculated from the primary loop. 

Secondary Regenerator and Bypass Valve Models 


The secondary regenerator is modeled as a two-port network. The 
empirical relationships used for this model are typical for a counterflow 
heat exchanger and were taken from Reference 7, The hot outlet tempera 
ture, t}i2t t>e written as 


^h2 = (^:^^h> "*^hl)'''thi (17) 

where 

Cc = Wc Cp(T) 

Ch = Wh Cp(T) 

W^, Wjj = coolant flow rates on cold and hot sides (Ib/hr) 

Cp(T) = coolant specific heat (a function of temperature) 

< = regenerator effectiveness 
t(;i =: cold side inlet temperature (“F) 
thi = hot side inlet temperature (“F) 
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Figure 4 is a diagram of the regenerator with assigned terms* 



AA A . 

* 



TTW * 

^ A A - - 

t 



' ■ 






Figure 4. Secondary Regenerator Diagram 


as 


From Equation (17), the capacity ratio, can be rewritten 



Wc AThAhc 


(18) 


where 


t- = cold-side outlet temperature (“F) 
^2 

AT = change in temperature ("F) 

Ah = change in enthaply (Btu/lb) 


Since the flow on the cold side of the regenerator is regulated by the 
secondary bypass valve, which is controlled by the primary condenser exit 
temperature, the portion of the cold-side flow that is bypassed must be 
considered. The relationship of the bypass flow rate, Wgp, and the bypassed 
fraction, a = Wgp/Wh, and the cold- and hot- side flow rates are given by the 

following equations: 


Wh = + Wgp 

(19) 

Wc/Wh == 1 - CK 

(20) 
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Equation (18) can then be rewritten by using the above expressions to give 


Cc/Ch = (1 -o) 


Ahh 


= P 


( 21 ) 


The parameters € and p, together with the factor to represent the 
time response, define the relationship between the inlet and outlet tempera- 
tures of the regenerator; 


tc2 


(1-Oni 

* ^12 


tci 

th2 


e (3X21 

(1 - 6p)x22 




( 22 ) 


This relation is solved in the program by a group of subroutines that 
update sliding arrays, move backward in the arrays to simulate delay, and 
perform integration. Figure 5 illustrates the relationship of the secondary 
regenerator cold-side inlet temperature to the regenerator effectiveness for 
various coolant flow rates. The data above an inlet of 70 degrees have been 
extrapolated. 

The secondary bypass valve characteristics are shown in Figure 6. 
The program keeps track of the TCE in ascending or descending order and 
correspondingly interpolates on the correct curve. 

Radiator Sygtem 


The fuel cell and radiator system consists of eight 5-sqoare-foot panels 
located on the CSM fairing. For the Skylab mission only two of the fuel cell 
and radiator coolant loops will be filled and connected for two- fuel-cell 
operation. While the radiator model discussion in this section is for the 
basic three-fuel-cell operational mode, two-fuel-cell operation for Skylab is 
simulated by using a zero flow pump characteristic for the third cell. For 
one-fuel-cell contingency operation, zero flow pump characteristics are used 
for two of the three cells. Reduced radiator area operation for low-power 
operation for Skylab missions is the same as for other Apollo missions; The 
bypass valve is actuated, and the last three panels are bypassed, and five- 
panel (five -eighths) radiator operation is achieved. 

The original radiator model was included with the data deck for the 
CINDA program received from NASA (Reference 8). The model included a 
nodal network for each of the eight radiator panels and the coolant delivery 
and bypass lines. Figure 7 illustrates the nodal model that is typical for a 
single radiator panel. The solid and fluid node and conductor numbers shown 
in this figure represent panel 1 in the radiator system model. The corre- 
sponding node capacitance and conductor valves are given in Table 4, The 
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GLYCOL TEMPERATURE = 130 F, PRESSURE 38 PSIG, TOTAL FLOW 80 PPH 



TEMPERATURE («F) 


Figure 6. Secondary Bypass Valve Characteristics 
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Table 4, Capacitance and Conductance Values 


Solid Node 

Solid Node Capacitance 

Solid 

Conductor Value 

Numbe r 

pCpV {Btu/“F) 

Conductor Number 

KA/L(Btu/Hr“F) 

701 

0. 043269 

701 

0. 199691 

702 

0, 037088 

702 

0. 176766 

703 

0, 043269 

703 

0, 090127 

704 

0. 098515 

704 

0. 121294 

705 

0. 080002 

705 

3. 111299 

706 

0.067608 

706 

3.927239 

707 

0. 067608 

707 

3.422037 

708 

0.080002 

708 

0. 106059 

709 

0. 067608 

709 

0. 182910 

710 

0, 067608 

710 

0. 072776 

711 

712 

0. 080002 
0.098515 

711 

0. 010989 



713 

0. 098515 

Radiation 

Conductor Value 

714 

715 

0. 080002 
0. 067608 

Conductor Number 

o-t FA(Btu/hr*F4) 



716 

0. 067608 

901 

0. 2298 X 10"9 

717 

0. 080002 

902 

0, 1970 X 10-9 

718 

0, 067608 

903 

0, 5234 X 10-9 

719 

0. 067608 

904 

0. 2700 X 10-9 

720 

0. 080002 

905 

0. 3591 X 10"9 

721 

0. 098515 



722 

0. 043269 



723 

0. 037088 



724 

0.043269 



893 

0.000000 




original values for the solid conductors were increased ten percent after it 
was found that the foil was two to three times thicker than required by the 
original specifications (Reference 9). Loss of heat through the radiator panel 
edges to the SM structure was accounted for by adding another solid node to 
the nodal network to represent the SM structure. Each of the solid-edge 
nodes of the eight radiator panels was tied to this SM structure node by a 
solid conductor. A structure temperature within a ±150 F range and an 
infinite capacitance is assigned the structure node, depending upon the 
environment being simulated. 
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The fuel cell /radiator coolant is a water- ethylene glycol solution 
composed of 62, 5 percent ethylene glycol by weight. The coolant properties 
used in the analysis as functions of temperature are density (Figure 8); 
conductivity (Figure 9); specific heat (Figure 10); viscosity (Figure 11); and 
relative enthalpy (Figure 12). The coolant pump characteristic, flow versus 
pressure drop in the radiator loop, is shown as Figure 13. The pressure 
drop in the radiator loop is calculated by using the Fanning equation with a 
dynamic head loss factor in the subroutine (Reference 8). 

Reactant Preheaters 


In the secondary loop, the water -glycol, after exiting from the con- 
denser, flows through the reactant preheaters before passing into the 
secondary regenerator. The basic equations simulating the reactant pre- 
heaters were taken from Reference 10, with corrections made to preserve 
compatibility with the secondary loop model. 

The effectivity of the oxygen preheater is shown in Figure 14 and 
expressed as 


EFF = 0. 53 + 0, 002 WDWG + 0. 018 ^ WD02 (23) 


where 


EFF = effectivity 

WDWG = water -glycol flow rate (Ib/hr) 

WD02 = oxygen flow rate (Ib/hr)- 

The oxygen preheater outlet temperature can then be found from 
Equations (23) and (24); 

T020 = T02I + (TWGOI-T02I) EFF (24) 

where 

T020= oxygen outlet temperature (^R) 

T02I = oxygen inlet temperature assumed constant := 530*R 
TWGOI = water-glycol temperature at condenser outlet (*R) 


- 25 - 


SD 70-266 






300 



SD 70-266 


Figure 12. Relative Enthcilpy o£ Water- Glycol Mixture 


space Division 

North American Rockwell 



OXYGEN FLOW RATE (LB/HR) 

(FROM P&WAMONTHLY PROGRESS REPORT, PWA - 2345, 25 MAY 1964) 

Figure 14, Oxygen Preheater Effectiveness 

From Equation (25), the change in enthalpy of the water-glycol flowing 
through the oxygen preheater can be obtained: 

DH02 = (T02I-T020)-(0. 152+0. 0000225*(T020+T021))>:< WD02/ WDWG (25) 

where 

DH02 = enthalpy change of the water-glycol through the oxygen 
preheater (Btu/lb) 

The temperature of the water -glycol (TWGHI) exiting the oxygen pre- 
heater is found from the change in enthalpy and preheater water -glycol inlet 
temperature. 

The effectivity of the hydrogen preheater is assumed constant at 
0.75 (Figure 15). Equation (26) is similar to Equation (24) and gives the 
hydrogen preheater outlet temperature; 

TH20 = TH2I + 0. 75 * (TWGHI-TH2I) (26) 
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(FROM P&WA MONTHLY PROGRESS REPORT, PWA - 2332, 27 APRIL 1964) 

Figure 15. Hydrogen Preheater Effectiveness 

where 

TH20 = hydrogen outlet temperature ( ®R) 

TH2I.= hydrogen inlet temperature, assumed constant = SSO^R 

TWGHI = water-glycol temperature at hydrogen preheater inlet (*R) 

The enthalpy change of the water -glycol flowing through the hydrogen 
preheater is then 

DHH2 = (TH2I-TH20) * 3, 47 *1* WDH2/WDWG (27) 

where 

WDH2 = hydrogen flow rate (Ib/hr) 

DHH2 = enthalpy change in the water -glycol through the hydrogen 
preheater (Btu/lb) 

The temperature of the water -glycol exiting the hydrogen preheater can 
then be found from the enthalpy change and the water -glycol temperature 
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entering the preheater* The water-glycol exit temperature is taken as the 
inlet temperature of the hot side of the secondary regenerator. 

Fuel Cell Heat Loss to Stru cture 

The fuel cell heat loss to the SM structure is based upon information 
given in Reference 5. The data as shown in Reference 5 are plots of fuel 
cell heat loss as a function of stack temperature for three parametric struc- 
ture temperatures: 30> 80 and 130 F* These data were replaced with the 

following linear equation, which is used in the SINDA fuel cell model: 

QSM == 3.25 TS - 1,6 .^= TA - 882, (28) 


where 


QSM = fuel cell heat loss to structure (Btu/hr per cell) 

TS = fuel cell stack temperature ( *F) 

TA = ambient structure temperature C^F) 

The ambient temperature is assigned as an independent parameter, 
based upon space environmental conditions. A temperature of 30 F is 
selected for extreme cold conditions; 1 30 F for extreme hot conditions, and 
80 F for nominal conditions. There is no attempt to determine the direct 
heat conduction through the fuel cell cone mount nor to determine the heat 
radiation from the pressure jacket and accessory package components. The 
program does not attempt to determine the changes in the ambient structure 
temperature resulting from such heating. The equation is accurate to within 
±4 percent of the data from Reference 5. 

4.2 PROGRAM DEFINITION 

The CINDA-3G (Chrysler Improved Numerical Differencing Analyzer 
for Third-Generation Computers) computer program was developed by the 
Thermodynamics Section of the Aerospace Physics Branch of the Chrysler 
Corporation Space Division at the National Aeronautics and Space Administra- 
tion’s Michoud Assembly Facility. The CINDA-3G program, written to run 
on the Univac 1108, was converted by the Computer Science Branch of the 
Idaho Nuclear Corporation to run on the IBM 360/75 computer. The con- 
verted program, called SINDA-3G (Systems instead of Chrysler), provides 
a variety of methods for the solution of thermal analog models presented in 
a network format. The network representation is unique in that it has a 
one-to-one correspondence to both the physical model and the mathematical 
model. The program allows the models to be developed through use of 
combinations of FORTRAN statements, user-initiated subroutines, and the 
numerous subroutines contained within the program. These program sub- 
routines can be used for handling interrelated complex phenomena such as 
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sublimation; diffuse radiation within enclosures ; simultaneous, one -dimensional, 
incompressible, fluid flow, including valving and transport delay effects; and 
similar areas associated with heat transfer and fluid flow. 

Li the hands of a competent engineering analyst, the SINDA-3G program 
is a powerful tool for analyzing thermal systems, 

4. 3 MODEL AND PROGRAM INPUT 

Data Deck Setup 


A SINDA-3G program deck contains two main blocks; a data block and 
an operations block. Each block is subdivided into four blocks. The four 
data blocks are entitled the NODE DATA, CONDUCTOR DATA, CONSTANTS 
DATA, and ARRAY DATA. The four operations blocks are designated 
EXECUTION, VARIABLES 1, VARIABLES 2, and OUTPUT CALLS. 

Card columns 12 through 80 comprise the data field. The instruction 
field (operations blocks) consists of columns 12 through 72. The program 
processes the problem data into FORTRAN common data and converts 
instructions into FORTRAN source language. They are then passed on to the 
system FORTRAN compiler. Instruction cards containing an F in column 1 
are passed on exactly as received. Discussion of the operations blocks 
will follow discussion of the data blocks. 

Data input to the data blocks may be one or more integers, floating 
numbers (with or without the E exponent designation), or alphanumeric words 
of up to six characters each. The reading of a word or number continues 
until a comma is encountered. Then the next word or number is read. Words 
or numbers nriay not be broken between cards, and a new card is equivalent 
to starting with a comma. Therefore, no continuation designation is 
required. When sequential commas are encountered, the program places 
floating-point zero values between them. Reading continues until the terminal 
column is reached or a dollar sign is encountered. Comments for a data card 
can be placed after a dollar sign and are not processed by the program. 

The first card in each of the four data blocks must be started in 
column 8 with the mnemonic code BCD (binary coded decimal), then an 
integer (1 through 9) in column 12, and the name of the data block starting 
in column 13. Each data block must be terminated with a mnemonic END 
card. 
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The NODE DATA block contains the data describing all nodes in the 
network. Each set of node data is grouped similar to the following code; 

column 12 

N#, Ti, Ca 


where 


N# = integer node number 
Ti = initial node temperature 
Ca = node capacitance 

All nodes are numbered by the program sequentially (from I on) in the 
order received. The user input number is designated the actual node number. 
The program assigned number is termed the relative node number. 

Second in the group of data blocks is the CONDUCTOR DATA block. 

In the nodal network, nodes are joined together with conductors. The three 
types of conductors used in the SINDA-3G program are conduction (solid), 
convection, and radiation conductors. Conductors are input to the program 
with the code. 


column 12 

G#. NA, NB, Cn 


where 


G# = integer conductor number 
NA = one adjoining node number 
NB = the other adjoining node number 
Cn = conductance value 

If more than one conductor has the same constant value, they may 
share the same conductor number and value. This is accomplished by 
placing two or more pairs of integer-adjoining node numbers between the 
conductor number and value. 

The data for the CONSTANTS DATA block are always input as doublets, 
the constant name or number followed by its value. They are divided into 
two types, control constants and user constants, and may be intermingled 
within the block. User constants receive a number. Control constants have 
alphanumeric names. 
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The ARRAY DATA block is the last of the four data blocks. Input of 
data into this block is exceedingly simple. The array number is listed^ 
followed by the sequential listing of the data^ and terminated with an END 
(data END, not mnemonic). The interpolation and matrix subroutines of the 
program make extensive use of these arrays. The SPACE option in the 
ARRAY DATA block is an easy way for the user to specify a large number of 
locations, which are initialized by the preprocessor program as floating- 
point zeros. When this option is used, the array number is listed, followed 
by the word SPACE and the number of locations to be initialized, and 
terminated with an END, 

While the four data blocks provide the data for the program, the four 
operations blocks determine the program control. This they do through the 
use of various operations and instructions. The four operations blocks, 
EXECUTION, VARIABLES 1, VARIABLES 2, and OUTPUT CALLS, are 
preprocessed by the SINDA-3G program and passed on to the system 
FORTRAN compiler as four separate subroutines entitled EXECTN, VARBLl, 
VARBL2, and OUTCAL, respectively. Figure 16 illustrates the basic flow 
diagram for the solution of the network* 

When the FORTRAN compilation is successfully completed, control is 
passed to the EXECTN subroutine. It sequentially performs the operations 
in the same order as input by the user in the EXECUTION block. None of 
the operations specified in the other three blocks is performed unless it is 
called for either directly by name in the EXECUTION block or internally by 
some other called for subroutine. All operations and instructions listed in 
EXECUTION block and performed by the program are executed only once. 
Because of this feature, the EXECTN subroutine can be used to initialize 
constants and variables, fabricate new arrays, establish steady-state 
parameters, and perform other operations that are completed only once 
during the duration of the program. 

The operations in the VARIABLES 1 block can be considered pre- 
solution operations. These operations may include construction of 
temperature arrays, establishment of heating rates, calculation of heating 
sources, or other basic operations required for solution of the thermal 
nodal network. 

In the same respect, the VARIABLES 2 block operations may be 
thought of as post-solution operations, VARIABLES 2 allows the user to 
look at the recently solved network. Typical operations of the VARBL2 
subroutine may include integration of flow rates, corrections of empirical 
relationships to reflect thermal solution of the nodal network, updating of 
conductances to account for changes in node temperatures, etc. 


- 34 - 


SD 70-266 



Space Division 

North American Rockwell 



OPERATION DESCRIPTION 


1 I 

I VARBLl ~1 

\ VARBL2~| 
I OUTCAL "I 
I MTC 

r~irn 


CALCULATE TIME STEP 
VARIABLES 1 OPERATIONS 
SOLVE NETWORK ■ 
VARIABLES 2 OPERATIONS 
OUTPUT CALLS OPERATIONS 
MODIFY TIME CONTROL 
ERASE ITERATION 


CHECK 

REVERSE DIRECTION IF 

<!> 

BACKUP NONZERO 


RELAXATION CRITERIA 
NOT MET 

<S> 

TIME OR TEMP CHANGE 
TOO LARGE 

<3> 

BACKUP NONZERO 

<3> 

NOT TIME TO PRINT 


PROBLEM STOP TIME NOT 
REACHED 


Figure 16. Basic Flow Chart for Network Solution Subroutines 
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The operations in the OUTPUT CALLS block are performed on the 
output interval specified by the user in the program* Since the operations 
are performed only at the output interval, OUTPUT CALLS typically con- 
tains only instructions for outputting information. 

The aforementioned data and operations blocks constitute a SINDA-3G 
program data deck. The deck must be terminated with the following card; 

^column 8 12 

BCD 3END OF DATA 

The user has the option to use the subroutines contained in the 
SINDA-3G library or to write his own. When non-SINDA-3G subroutines are 
being called, the data communication is obtained through subroutine argu- 
ments similar to any other subroutine. 

Sample Problem Input Data 


The listing of the sample input data is included in the appendix. The 
data are presented in the same order as stated in the Data Deck Setup section; 
the node data first, followed by the conductor data and constants data, and 
concluded with the array data. 

The node data includes the solid temperature nodes for the eight 
radiator panels, the water -glycol fluid temperature nodes for the three- 
fuel-cell system, the fluid pressure nodes, and the edge sink and environ- 
mental sink nodes for the eight panels. 

The 11 solid (conduction) conductors are the first of the data listed in 
the conductor data. Next are the convection conductors and the pressure 
and fluid flow conductors for the three systems. The conductor data are 
concluded with the five radiation conductors. 

The constants data contain the 472 constants used in the program* 

These constants are listed by their actual number rather than the program- 
assigned relative numbers. 

Listed in the last section of the problem input data are the array data. 
Eighty-four arrays are listed, with most arrays containing data and the 
remaining arrays using the SPACE option to allocate program storage 
locations . 

Following the input data is a listing of the main program (EXECTN, 
VARBLl, VARBL2, and OUTCAL) and the 12 user-written subroutines* 
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Sample Problem Output 

Figure 17 illustrates a typical problem output listing. The listing is 
headed with the title "Systems Improved Numerical Differencing Analyzer ^ 
StNDA ^ North American Rockwell Corporation - Space Division" and 
followed by 11 output interval data sets. Each output interval data set is 
started with four asterisks. The first data line contains five control con- 
stants with their values listed to the right of each constant. The last three 
constants have the relative node number enclosed in parentheses. The 
control constants are; 

Present mission time (hr) 

Last time step used for transient network problem (hr) 

Most recent minimum stability criteria for the 
network 

DTMPCC Maximum diffusion temperature change calculated 

over the last time step 

ARLXCC Maximum arithmetic relaxation change calculated 

over the last iteration 

The next four lines of output are associated with the primary loop 
operation. The first two lines list the titles^ with the values following. 

The titles are defined as follows: 


TIME 

DTIMEU 

CSGMIN 


SYSTEM Fuel cell system in operation (for Skylab configuration 

only two fuel cell systems are used, systems 1 and 3) 

POWER Fuel cell power output (watts) 

CURRENT Fuel cell output current (amps) 

VOLTAGE Fuel cell output voltage (volts) 

TSI Stack inlet temperature (*F) 

TSE Stack exit temperature (®F) 

TCIP Condenser inlet temperature, primary loop ( *F) 

TCEP Condenser exit temperature, primary loop(“F) 

QCOND Heat flow across the condenser (Btu/hr) 
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W-COND Weight rate of water condensed (Ib/hr) 

WS-RATE Percent of water produced in the time interval that 
remains in the stack {%) 

Oxygen reactant temperature at stack inlet (“F) 

Hydrogen reactant temperature before mixing with the 
primary recirculation stream (“F) 

PCKOH Electrolyte concentration ratio 

Secondary loop operations are listed in the next four lines, with the 
titles followed by their values. The titles represent: 

FLOW RATE Water glycol flow rate (lb /hr) 

QRAD Heat flow across the radiators (Btu/hr) 

TRADIN Temperature of the radiator inlet {“F) 

TRADOUT Temperature of the radiator outlet { ®F) 

DPRAD Pressure drop through the radiator loop (psi) 

TSRCI Temperature of the secondary regenerator at the 

cold-side inlet (®F) 

TSRCE Temperature of the secondary regenerator at the 

cold-side exit (“F) 

TSRHI Temperature of the secondary regenerator at the 

hot-side inlet {"F) 

TSRHE Temperature of the secondary regenerator at the 

hot-side exit { “F) 

TCIS Condenser inlet temperature, secondary loop ( "F) 

TCES Condenser exit temperature, secondary loop (°F) 

BPFS Bypassed ratio at the secondary regenerator 

The final data group lists the radiator inlet and outlet temperature for 
the eight radiator panels for all three systems. The temperatures are 
listed in degrees Fahrenheit. The four asterisks following this data group 
designate the start of a new output interval data set. 
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4,4 ANALYTICAL RESULTS 

Three temperature levels are monitored to determine if the fuel cells 
are operating within the specified limits. These temperatures are measured 
at the condenser primary exit^ the radiator exit, and the stack. The nominal 
condenser primary exit temperature range is 155 to 165 F, with allowable 
cycling to 200 F, The caution and warning (C&W) alarms for the condenser 
exit are set for 150 F and 175 F, For the analysis, the acceptable minimum 
power operation was defined for condenser operation as between 150 and 
155 F. The acceptable maximum power operation was defined as between 
165 and 200 F* The stack and radiator exit temperatures are less limiting 
than the condenser exit temperature. This is evidenced by results in which 
seven of the eight cases reached the condenser exit limits. The remaining 
case was limited by the stack temperature. 

The nominal stack temperature range is 390 to 460 F, with the C&W 
alarms at 360 F and 475 F, For the analysis, the acceptable minimum stack- 
temperature range was 360 to 390 F, The maximum was defined as 460 to 
475 F, The radiator exit-temperature nominal range is 0 to 120 F, with 
allowable high-temperature cycling to 180 F, The C&W alarm is set only 
for the minimum temperature condition of -30 F, Only the acceptable 
minimum radiator temperature range was defined for the analysis. This 
range, -30 to 0 F, was not encountered in any of the eight cases. 

Two environmental heating profiles were used for the heat that is 
radiated to radiator panels. The hotter environment, which was used for 
the maximum power cases, has a p angle of 73, 5 degrees for a circular 
earth orbit of 235 nautical miles with the vehicle in a Z -local vertical 
attitude hold. The colder environment for the minimum power cases uses 
a p angle of 0 degrees for the same orbit and attitude conditions. 

The power levels shown in Tables 1, 5 and 6 represent the expected 
power requirements for the Skylab missions. One-fuel-cell operation and 
five-eighths radiator area are contingency operation modes that are 
considered irreversible. That is, once the contingency mode is selected, 
a return to normal two-fuel-cell or full -area operation is never required. 

With the exception of Case 2, the power levels shown in Table 5 can 
be raised slightly by operation at the condenser C&W limit of 175 F, 
Operation to 200 F is acceptable for peak cycling conditions. However, once 
the nominal bypass limit of 165 F is exceeded, a slight increase in power 
results in a significant temperature rise. Case 3, one fuel cell, full area 
at 60 amperes, is the only case where the stack temperature limits the power. 
In this case, a greater increase in power is available if operated to the stack 
temperature C&W limit of 475 F, 
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For all of the minimum power cases listed in Table 6, slightly lower 
power capability is possible by operating down to the condenser C&W limit 
of 150 F, The temperatures shown in Tables 5 and 6 are orbital extremes 
except where differences occurred between fuel cells in two-fuel-cell opera- 
tion. Then the value shown is the average between the two cells at the 
orbital extreme condition. 


Table 5. System Temperatures for High Power Levels 


Case 

p 

Angle 

(deg) 

Number of 
Fuel Cells in 
Operation 

Radiator 
Area {operating 
panels) 

Total Fuel Cell 
Current 
Level (amps ) 

System T em 

peratures 

rFi 


Stack 

Radiator 

Inlet 

Radiator 

Outlet 

n 

BB 

2 

8 

100. 

168. 45 

453. 93 

189. 87 

150. 58 

mm 

wSB 

2 

5 

90. 

192. 01 

452. 11 

208. 98 

178. 48 

B 

mB 

1 

8 

60. 

163. 77 

470. 24 

190. 24 

L38. 04 

B 

B 

1 

5 

50. 

172. 98 

458. 70 

194. 06 

[55. ^ 1 


Table 6. System Temperatures for Low Power Levels 


Case 

P 

Angle 

{degl 

Number of 
Fuel Cells in 
Operation 

Radiator 
Area (operating 
panels) 

Total Fuel Cell 
Current 
Level (ampsl 

Syst 

em Temperatures 

( ’F) 

Condenser 

Exit 

Stack 

Radiator 

Inlet 

Radiator 

Outlet 

5 

0. 

2 

8 

50. 

154.30 i 

418. 25 

52. SB 

33. 03 

6 

0, 

2 

5 

40. 

151.46 ! 

412. 36 

66. 1 1 

40 . 82 

7 

0. 

1 1 

8 

30. 

155.06 ' 

426. 05 

47. 12 

19. 08 

8 

0. 

1 ' 

5 

25. 

153.41 

419. 54 

46, 71 

24, 57 
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APPENDIX: INPUT DATA, MAIN PROGRAM LISTING, 
AND USER SUBROUTINE LISTINGS 


INPUT DATA preceding PAGE BLANK NOT FILMED 

SCO SPCS 

FNri 

SCO 3F40DF DATA 

RFM radiator solid tempfrature nodes IPANEL U 

7 31, 103. t .0432693. 702, IOC. . .037088C.703 , 100, , .C4 326Q3 

734.100. , .0985149,705, 100, , . 08000 1 8, 7C6 , 1 00 . , .067608 3 

737.100.. .0676083. 708. 100... 0800018. 709. 100. . .067608} 

71 0.100.. . 06 76083.711. 100-,, 080001 8, 71 2, 100. ,.0985149 

713.103.. .0985149. 714. 100- , .030001 8 , 71 5 , 100. , .06^6083 

716.100.. .067608 3. 7 17. 1 00... 08000 1 8. 71 8. 100. . .0676083 

719.100.. .0676083. 720. 100... 08000 18. 72 1. 1 00. ..0985149 

722.100., .0432693, 723, 100-,. 0370880, 724, 1 0.,. 0 32. ' 

RFM RADIATOR SOLID TEHPFRATUP6 NODES (PANEL 21 

725.100.. .0432693. 726. 1CO.,. 03 70880, 727, IOC. , .0432693 

728. 100.. . 0985149. 729. 100... (5800018, 730, 1 00.,. 0676C83 

731. 100.. . 0676083. 732. 100... 0800018.733, toe. ,.0676083 

734.100.. .0676083.735, too.,. 08000 18, 736, 100. ,.0985149 

737. 100. . .0985149.738. 100. ..080001 8. 739. 100. . .0676083 

740. 1 00.. . 0676083. 741. 1 00... 080001 8. 742. 100. ..0676083 
743, too... 0676083, 744, ICO*,. 08000 18, 745, 100. ,.0985149 

746. 100.. . 0432693. 747. 100... 0370880. 748. 100. . .0432693 
RFH RADIATOR SOLID TEMPERATURE NODES (PANEL 31 

749. 100.. .0432693.750. 100*,. 0370880, 751 , 100. , .0432693 
752,100*, .09 851 >,9, 753 , 100. , .O.IOC 1 8 , 754 . 1 DC . , .0 ^ 7 60' 1 
755,100*, .0676083,756, too* , .0800018,757, 100. ,.0676083 

758. 100. , .0676083,759, 100.,. 080001 8, 760,100. , .0985149 

761. 100.. . 0985149. 762. 100... 0800018.763. 100... 0676083 

764.100. . .0676083.765. 100. . .080001 8 . 766 . 1 00. . .0676083 

767. 100.. . 0676083. 768.100... 08000 18. 769. 100. . .0985149 

770.100.. .0432693. 771. 100... 0370880. 772. 10C. , .0432693 
REM RADIATOR SOLID TEMPERATURE NODES (PANEL 41 

773. 100. . .0432693. 774. 100. ♦ .0370880, 775 , 100. , .0432693 

776. 100.. . 0985149. 777. 100... 080001 8. 778, too. ,.0676083 
779, 100*, *06 76083, 780, 100.,. 080001 9, 78 1, 100.,. 0676083 

792. 100.. . 0676083. 783. 100... 080001 8. 784. 100... 0985149 

795.100.. . 0985149. 786. 100... 080001 8. 767. 100. ..0676083 

799.100.. .0676083. 789. 100... C80001 8, 790, 100. ,.0676083 
791,lC0«,ii0676083,792, 100.,.C800018,T93, 100. . .09 85149 

794. 100. . .0432693. 795. 100.. .0370880. 796 . 1 00. . .0432693 

REM RADIJtTOR SOLJO TEMPERATURE MOOES (PANEL SI 

797^100. ',.0432693, 798,100.,. 0370880, 799, too.,. 0432693 
890, too., *0985149, 801, 100.,. 080001 a, 802, 100.,. 0676083 

893. 109., *0676083, 004, 100. ,.0800018,805,100. , .0676083 

806. 109.. . 0676083. 807. 100... 080001 6. 805. 1 00. . .0985149 

809,109. 'I: .f.%^ L8,Dll,i. ‘ 

812.100., .0676083,8! ’.,100. ,.080001 8, 81 4, mo. , . > ■'AO" 

81 5. 109.. . 0676083. 816, too... 080001 8, 8 1 7, 100. ,.09851 49 

818.100.. .0432693. 819. 100*, .03 70880, 820, 100. , .0432693 
REN RADIATOR SOLID TEMPERATURE MOOES (PANEL 61 

821.100.. .0432693. 822. 100. ..03 70880. 823. 100.. .0432693 
824, too.,. 0985149, 825, too.,. 0300018, 826, 100. ,.0676083 

827. 100. *.0676083, 828, 190.,. 08000 18, 829, IOC.,. 0676083 

830. 109.. . 0676083. 831 .100...08000 18. 832. 100... 0985149 

833.109. . .0985149. 834. 190.. .080001 8. 835.100. ..0676083 

836. 109.. . 0676083. 837. 100... 089001 8. .838. 100... 0676083 

839.100.. .0676083.840.100. ..0800018.841.100. ..0985149 

842. 100 .. .0432693. 843. too., . 0370880. 844. 100. , .0432693 

Preceiiiig page Uank . 4, . 
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irrLOOiD 

tMDDOOlC 

$NOD0020 

tNODODTC 

*Nno0040 

tNODOOSC 

$NOD0060 

tNODOOTO 

tNODOCPC 

INODOC9D 

ANODOICO 

4NODOUO 

tNOOOl20 

tNODCUu 

tNOOOlAO 

tNDOOlSO 

$NDD0160 

SN0D0170 

SNOD018C 

tNDD0190 

TN000200 

tNOD02lO 

tN-100220 

INa 00230 

SNOOC240 

ANOD0250 

SNOD0260 

tNOn 027 O 

»NQD0280 

$NQD0290 

TNODOBOO 

SN000310 

SNOD0320 

tNOD0330 

$NOD0340 

ANOD0350 

SN000360 

SNODO370 

«N000380 

ANOD0390 

8N000400 

$NOD0410 

SNODC420 

<MU00430 

$NOD0440 

SN000450 

SN000460 

tN000470 

»M0D048.0 

$NOD0490 

SNOOOSOO 

SNOD0510 

$N000520 

SN000530 

SNOD0540 

$N000550 
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(IGM BAOIATOft SOLID TEMPERATURE NODES ("ANEL 7t 

S'tSt 103. » .0432h93,84f», 1 00 ... 0370880. 847, 100. , ,0432693 
848. ICO. . .0965149.849, 100... 083001 8.650.100. *.0676083 

8 51. 100.. . 0676083. 8 52. 1 00... 080 00 1 8. 853.1 00... C67608 3 

654.100.. .0676083.855. 1 00. . . 080001 6 . 856 . 1 00. . . 09 85 1 49 

8 5 7. 100.. . 0985149.858. 100... 080001 8. 859. 100. ..0676083 

360. 103.. . 067608 3. 86 1.130... 080001 8. 862. 1CC... 0676083 
86 3. 103. , .0676083.864. 100. .. 08000 1 8 . 865 . 100. ,. 09 851 49 

866.100. . .0432693.867. 100. .. 0370880. 868 . 100. .. 0432693 
REM RADIATOR SOLID TEMPERATURE NODES (PANEL 8) 

869.100.. .0432693.870* 100 ... 0370880. 871 , 100. . .0432693 
872. too... 09851 49, 873,100*,. 0800018. 874, 100... 0676083 

875. 100.. . 0676083. 876. 1CO.,. 0800018.677. 100... 0676083 

878.100. . .0676083. 879.100. ..0800018.880, LOO. ,.0985149 

881.100.. .0985149.882.100. ..0600018.883.100. . .0676083 

884.100.. .0676083.885.100. ..0800018.886. 100. ..0676083 

897.100. . .0676083.888.100. ..0800018.889. 100. ..0985149 

890.100. . .0432693. 891. 100.. .0370880.892.100. . .0432693 
REM FLUID TEMPERATURE ARITHMETIC NODES 

617.116.. -l.,616,119..-1..6l5.119..-l.,614,122..-l. 
6l3,t22.,-l.,612.L2S.*>l..611.12S.,-l.,610,128.,>l. 

609.128. . -l. .608. 131. .-1. . 607* 1 3 1 . , -1 . . 606 , 1 34 . 1 , 

605. 1 34.. -1.. 604. 137., -I,, 603. 1 37., -I,, 6C2, 140.. -1. 

635.116.. -l.,634,L19..-l.,633,119..-1..632,l22..-l. 

631.122.. >I..630.l25..-l.,629.l25..-l.,62e,L2A.,-l. 
62T,129.,-L..6Z6,13t.,-l.,625,13l.,-1..624,134.,-l. 
623.134.fl.,622,l3T..-1..62l,137..-1..620.140..-l. 

653.116.. -l.,652.U9..>l.,65l,119,.-1..650.122..-l. 

649. 1 2 2.. -1., 648, 125.. -1., 647, 1 25.. -1., 646, 12 8., -I . 

645.128.. >1..644.131..-1..643,131,.-l.,642,l34. ,>1. 

641.134.. >1.,640.137.,-1.,639,137..-1.,638.140.,-1. 

RFM FLUID TEMPERATURE NODES - M1RINS (SYSTEMS 1.2.3) 

6 18. 1 1 6.. -1., 636, 116., -1.. 654. 11 6., -1. 

REM FLUID TEMPERATURE- NODES - BOUNDARY (SYSTEMS 1,2.3) 

•601 ,l90...-6l9.190..,-637, 190., 

REM FLUID PRESSURE NODES (SYSTEM 1) 

-S0l,.,-502.,.-503,,,-504.,,-505..,-506,.,-507.,,-508.. 
-509.,,-510.,,-511.,,-512,,,-5U,..-5l4,,,-515,..-516.. 
•517. ..•518, . .•519. ..-520.. . “52 1.. .-522. ..-523* . .-524. . 
-525, ,.-526, ,.-527,0. 

REM FLUID PRESSURE NODES (SYSTEM 2) 

-528... -529... -530... -531.,. -532... -533... -534... -535,. 
-536. ..-537. t .-538. . .•539. . . -540* . .“54 1...-542. . . -543 . . 
-544.,, -545... -546... -547,., -548... -549... -550... •551,, 
•552. ..—553. . .—554.0. 

REM FLUID PRESSURE NODES (SYSTEM 3) 

-555... -556.,. -557... -SSS... -559... -560,,. -561... -562., 

-553 . . . -554. ..—565. . .—566. ..-567.. .-568.. .-569. . .-570. . 
-57l,;,-572,,.-573.,,-574.,,-575,,,-576...-577.,,-578.. 
-579,,,-5B0,, ,-581,0. 

REM EDGE SINK TEMPERATURE NODE FOR PANELS 1 TO 8 
-893.-150.. 

REM ENVIRONMENTAL SINK TEMPERATURE NODES FQR PANELS 1 TO B 

-901, -460... -902,-460.,, -903, -460.,, -904, -460... -905, -460.. 
-906, -46().,. -907, -460.., -908. -460., 

FNO 

BCD 3CQN0UCT0R DATA 

REM 

REM 

REM SOLID CONDUCTOR - NUMBER 701 
REM 

701,701,702,702,703,722,723,723,724 SPANEL I 

725,726,726,727,746,747,747,748 5PANEL 2 

749,750,750,751,770,771,771,772 SPANEL 3 


ANOD0560 

tN000570 

*NODO5B0 

iNOD0540 

tNOnO60O 

sNuno5io 

tNODO520 

iN0D0630 

tN0DO640 

tNOD0550 

SNODOBSC 

tNO0C670 

SN0D0680 

*NOD0590 

SN000700 

$NnO07l0 

$NODO720 

AN000730 

$N00075C 


tN0D0880 

SNOOD900 

SN0O09L0 

$N000920 

«N000930 

SNOD0940 

tN0D0950 

$NOD0960 

SN0D0970 

kNOD0980 

INOD0990 

SNOOIOOO 

tNODlOlO 

$N0Dt020 

SN001030 

tN0D1040 

SNODIOSO 

$N0D1060 


SN001070 

SNODIOSO 

SNOD1090 

tNODlllO 

ICNDOOIQ 

SCND0020 

tCND0030 

$CN00040 

$CN00050 

$CN00060 

$CN00070 

$CN00080 
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773, 77<t, 77^,,7 75,794, 7<J5, 795, 796 
79 7,798,79(3,799,813,819,819,820 
821,822,822, 823, 862, 863, 843,846 
865,866,866,867,866,867,867,668 

869.870.870.871.890.891.891.892, 
RFM SOLID COMOUCTOR - NUHBFR 702 

702,701,706,703,712,713,722,721 

725.726.727.736.737.766.765.768 

769. 752.751.760.761.770.769.772 

773.776.775.786.785.796.793.796 
797,800,799,808,809,818,817,820 

821.826.823.832.833.862.861.866 

865.868. 867.856.857.866.865.868 

869.872.871. 880.801 .890.889.892, 
RFM SOLID CONDUCTOR - NUMBFR 703 

703,701,705,702,700,703,711,716 
725,729, 726, 732,727,735, 738,766 
769,753,750,756,751,759,762,770 

773.777.776.780.775.783. 786.796 
797,801,798,806,799,807,810,818 
821,825,822,828,823,831,836, 362 

865.869.066.852.867.855.958.866 
869,873,870,876,871, 979,882,890, 

RFM SOLID CONOUCTDR - NUMBER 706 

706, 701 , 706,702 , 707, 702 , 709, 703 
715,722, 716,723,718,723, 719,726 
725,730, 726,731,726, 733, 727, 736 
7 39 , 766 , 760 , 76 T , 76 2 , 767 , 763 , 76 8 

769.756.750.755.750.757.751.759 

763.770.766.771.766.771. 767.772 

773.778.776.779.776.781.775.782 

787. 796.708.795.790.795. 791. 796 
797,802,799,803,798,905,799,006 
811,818,812,819,-816,819,815,820 

821 .826. 822.827.822.829.823.830 

835.862.836.863.830.863.839.866 
865,850,866,851,866,853,867,856 

859.866. 860.867.662.867.863.868 
869,876,870,875,870,877,871,878 

883.890.886.891.886. 891.887.892, 
PFM SOLID CONDUCTOR - NUMBER 705 

705.706.705.711.712.713.716.720 

728.729.735.736.737.738.766.765 

752.753.759.760.761.762. 768.769 
776,777,703,796,785,736, 792,793 
803,801, 907,808,809,810, 816,817 
826,825, 831, 832,833,836, 860,841 

868.869.855.856.857.858.866.865 
872,873,879,380,881,882,888,889,3 

RFM SOLID CONDUCTOR - NUMBER 706 

706,705,706,707,708,708,709,710, 

716.715.716.717.717.718.719.720 
729,730,731,732,732,733,736,735 

738.739.760.761.761.762.763.766 

753.756.755.756.756.757.758.759 
762,763, 766,765,765,766, 767,768 

777.778.779.780.780.781.782.783 
786,787, 788,789,789,790, 791,792 
801 , 802,803,806,806, 805. 806,807 
810,811,812,813,813.816,815,816 

925.826.827.828.828.829.830.831 
836,835,036,837,837,038,839,860 
869,850,851,852,852,853,856.655 

858.859.860.861.861.862.863.866 



3P6NFL 6 

SCNDOnqc 


iPANEL 5 

SCNDCIOC 


SPANEL 6 

SCNOOllC 


3PANEL 7 

SCND0120 

. l«9691 1 

$PANEL 8 

SCN00130 

SCN00160 

,726 

SPANEL 1 

tCNDOlSO 


tPANEL 2 

tCNDC-l6C 


$PANEL 3 

SCND017C 


SPANEL 6 

SCNDOIRO 


SPANEL 5 

SCNDOi90 


SPANEL 6 

1CND0200 


SPANEL 7 

SCN00210 

.1767653 

IPANEL 8 

SCNOG220 

ICNDO230 

, 722,717, 

723,720,726 

SPNL 1 SCNDC26C 

,761,767, 

766,768 

iPNL2SCNl>0?50 

,765,771 , 

768,772 

SPNL3SCND0260 

,789,795, 

792,796 

SPNL6SCND0270 

,813,819, 

816,020 

SPNL5SCND0280 

,837,863, 

860,866 

SPNL6SCN00290 

t S6 1 1 ^ 6 7 f 

866,868 

SPNL7SCND0300 

885, 891, 888, 892,. 090 127 

SPNL8tCND0310 

tCN00320 

*710 $ 

PANFL 1 

SCND0330 


PANEL 1 

SCND0360 


PANEL 2 

SCN00350 


PANEL 2 

SCND0360 


PANFL 3 

SCN00370 


PANEL 3 

SCND0380 


PANEL 6 

SCN00390 


PANEL 6 

tCND060C 


PANEL 5 

SCND0610 


PANEL 5 

SCN00620 


PANEL 6 

SCND0630 


PANEL 5 

SCND0660 


PANEL 7 

SCN00650 


PANFL 7 

SCND0660 


PANEL 0 

SCND0670 

.121294 $ 

PANEL 8 

SCND0680 

tCND0690 

*721 

SPANEL 1 

SCN00500 


SPANEL 2 

SCND0510 


SPANEL 3 

SCN0052C 


SPANEL 6 

SCND0530 


SPANEL 5 

1CND0560 


SPANEL 5 

SCND0550 


SPANEL 7 

SCND0560 

111299 

SPANEL 9 

$CN00570 

SCNOOSBO 

fTll 

SPANEL 1 

SCN00590 


SPANEL 1 

SCND060Q 


SPANEL 2 

tCN006lO 


tPANEL 2 

SCN00620 


SPANEL 3 

SCN00630 


SPANEL 3 

SCN00660 


SPANEL 6 

SCN00650 


SPANEL 6 

1CND0650 


SPANEL 5 

SCN00570 


SPANEL 5 

SCN00580 


SPANEL 6 

SCND0590 


SPANEL 6 

SCND0700 


SPANEL 7 

SCN00710 


SPANEL 7 

SCN00720 
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37 3,«7<i,«75*87fi»876, 877, 878,87<> 

88?, 883, 834, 885, 885, 886, 887, 888, 3. Q27239 
SOHO CONNECTOR - NUMBER 707 
70 7,706,707,709,7X0,715, 716,718,7X9 

730.731.733.734.739.740.742.743 

754.755. 757.758.763.764. 766.767 

778.779.781.782.787.783.790.791 
802,803,805,806,811,812,314,815 
826,827, 829,830, 835,836,838,839 
850,351,853,854,859,860,862,863 

874,875,877,878,883,884,886,887,3.422037 
REM SOLin CONDUCTOR - NUMBER 70S 

708,704,713,712,721,728,737,736,745 

752,761,760,769,776,785,784,793 

800,809,808,817,824,833,832,841 

843. 857. 856. 865. 872. 88 1.880. 889.. 106059 
RFM SOLID CONDUCTOR - NUMBER 709 

709,705,714,708,717,711,720 

729,738,732,741,735,744 

753,762,756,765,759,768 

777, 786,780,789,783,792 

801,810,804,813,807,816 

825,834,828,837,331,840 

849,858,852,861,855,864 

873.882. 876.885.879.888.. 1829098 iPANEL 
REM SOLID CONDUCTOR - NUMBER 710 

713,706,715,707,716,709,718,710,719 

730.739.731.740.733.742.734.743 

754.763.755.764.757.766.758.767 

778.787.779.788.781.790.782.791 
802, Sir, 803, 612, 80 5, 814, 806,815 
826,835, 827,836,829,838, 830*839 
850 , 859, 851 , 860,853*862, 854, 863 

874.883.875. 884. 877.886.876.887. . 0727762 
REM SOLID CONDUCTOR - NUMBER 711 

711,701,893,702,893,703,893,712,893,721, 

722,893,713,893,704,893 


(PANEL 

8 

(CND0730 

(PANEL 

8 

(CND0740 

(CND0750 

(PANEL 

1 

(CND0760 

(PANEL 

2 

(CND0770 

(PANEL 

3 

(CN00780 

(PANEL 


tCNn0790 

(PANEL 

5 

sCNDoeco 

(PANEL 

6 

(CND0610 

IPANEL 

7 

(CN00820 

(PANEL 

8 

(CN00830 

ICND084C 

(PANELS 

1-2 

(CND0R50 

(PANELS 

3*^ 

(CNDO06O 

(PANELS 

5-6 

(CND0870 

(PANEL 

8 

(CNDOBRO 

(CND089C 

iPANEL 

1 

(CND0900 

(PANEL 

2 

(CNDC910 

(PANEL 

3 

iCN0092C 

(PANEL 

4 

(CN00910 

(PANEL 

5 

(CND0940 

(PANEL 

6 

(CN00950 

(PANEL 

7 

(CND0960 

(CND0970 

(CN00980 

(PANEL 

1 

(CND0990 

(PANEL 

2 

(CNOIOOC 

(PANEL 

3 

(CNDIOIC 

(PANEL 

4 

(CN01020 

(PANEL 

5 

(CND103Q 

(PANEL 

6 

(CN0104C 

(PANEL 

7 

(CND1050 

(PANEL 

d 

ICN01060 


893,724,893,723,893 
(PANEL 1 


725.893.726.893.727.893.736.893.745.893.748.893.747.893 

746.893.737.893.728.893 t PANEL 2 

749.893.750.893.751.893.760.893.769.893.772.893.771.893 

770.893.761.893.752.893 » PANEL 3 

773.893.774.893.775.893.784.893.793.893.796.893.795.893 

794.893.785.893.776.893 t PANEL 4 

797, 893, 796,' 893, 799, 893, 808, 893, 8 17, 693, 820, 893, 819, 893 

818.893.809.893.800.893 t PANEL 5 

821.893.822.893.823.893.832.893.841.893.844.893.843.893 

842.893.833.893.824.893 S PANEL 6 

845.893.846.893.847.893.856.893.865.893.868.893.867.893 

866.893.857.893.848.893 $ PANEL 7 

869.893.870.893.871.893.880.893.889.893.892.893.891.893 

890. 893. 881. 893. 872. 893. . 0109890 ( PANEL 8 

REH 

REM CONVECTION CONDUCTORS 
REM 

1001. 602. 705.. 1002. 603. 71 7.. 1003 .604. 735.. 1004. 605. 744, 

1005. 606. 753.. 1006. 607. 762.. 1007. 608. 780.. 1008. 609. 789, 
1009* 610, 807,, 1010, 61 1,816,, 10 11, 6 12, 825,, 1012, 61 3, 834, 

1013. 614. 852.. 1014. 615. 861.. 1015. 616. 879.. 1016. 617. 866, 

1017. 620. 71 1.. 1018. 621. 720.. 1019. 622. 729.. 1020. 623. 738, 

1021. 624. 756.. 1022. 625. 765.. 1023. 626. 783.. 1024. 627. 792, 

1025. 628. 801.. 1026. 629. 810.. 1027. 630. 828.. 1028.631.837, 

1029. 632. 855.. 1030. 633. 864. .1031 .634. 873. . 1032.635.882 , 

1033. 638. 705.. 1034. 639. 714.. 1035. 640. 732.. 1036.641. 741, 

1037. 642. 759.. 1 038. 643. 768.. 1039. 644. 777.. 1040. 645. 786, 


(CN01070 

3CN01080 

3CND1090 

tCNDllOO 

ICNDlllO 

$CN01120 

ICNOinO 

tCNOlUO 

SCN01150 

$CN0U60 

(CNDllTO 

tCNOliaO 

$CND1190 
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North American Rockwell 


PEM 

REM 

REM 


REM 

REM 

REM 


REM 


804 1 , 1042,647, 813 1043,648, R31 , 1 10^^*649,840* 
1045, 650, 849,, 1046, 651, 858,, 1047, 652, 876,, 1048, 653, 885, 

PRESSURE CONDUCTORS (SYSTEM 11 

501.501.502, 

502.502.503, 

503.503.504, 

504.504.505, 

505.505.506, 

506.506.507, 

507.507.508, 

500.508.509, 

509.509.510, 

510.510.511, 

511.511.512, 

512.512.513, 

513.513.514, 

514.514.515, 

515.515.516, 

516.516.517, 

517,517,510, 

518.518.519, 

519.519.520, 

520.520.521, 

521.521.522, 

522.522.523, 

523.523.524, 

524.524.525, 

525.525.526, 

526.526.527, 

527,516,526, 

PRESSURE CONDUCTORS (SYSTEM 2) 

529.528.529, 

529.529.530, 

530.530.531, 

531.531.532, 

532.532.533, 

533.533.534, 

534.534.535, 

535.535.536, 

536.536.537, 

537.537.538, 

538.538.539, 

539.539.540, 

540.540.541, 

541.541.542, 

542.542.543, 

543.543.544, 

544.544.545, 

545.545.546, 

546.546.547, 

547.547.548 , 

548.548.549, 

549.549.550, 

550.550.551, 

551.551.552, 

552.552.553, 

553.553.554, 

554,543,553, 


tCNOUOC 

$CND1210 

tCND12?0 

*CND123C 

tCN(U24C 


tCN0126C 

$CND1270 

tCN01280 

$CND1290 


SCN01310 

SCN01320 
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PfM PRESSUftf: CONOLtCTflRS (SVSTEM 3) 
REM 

555,555*556, 

556.556.557, 

557.557.558, 

558,558,556, 

556.556.560, 

560.560.561, 

561.561.562, 

562.562.563, 

563,563,566, 

564.566.565, 

565.565.566, 

566.566.567, 

567.567.568, 

568,568,566, 

569.568.570, 

570.570.571, 

571.571.572, 

572.572.573, 

573.573.574, 

574.574.575, 

575.575.576, 

576.576.577, 

577.577.578, 

578.578.579, 

579.579.580, 

580.580.581, 

581,570,580, 

REM 

REM FLUID FlOH CONDUCTORS (S7STEM II 
REM 


601 

,-601 

,602, 

602 

,-602 

,603, 

603 

,-603 

,604, 

604 

,-604 

*605, 

605 

,-605 

,606, 

606 

,-606 

, 607 , 

607 

,-607 

,608, 

60S 

,-608 

, 609, 

609 

,-609 

,610, 

610 

,-610 

,611, 

6U 

,-611 

,612, 

612 

.-612 

,613, 

613 

,-613 

,614, 

614 

,-614 

,615, 

615 

,-615 

,616. 

616 

,-616 

,617, 

617 

,-617 

,618, 

618 

,-611 

*618. 


REM 

REM FLUID FLOW CONDUCTORS (SYSTEM 21 
REM 

619, -619,620, 

620, -620,621, 

621,-621,622, 

622, -622,623, 

623, -623,624, 

624, -624,625, 

625, -625,626, 

626, -626,627, 

627, -627,628, 

628, -628,629, 

629, -629,630, 


eCN(U330 

tCNfH340 


$CN01360 

$CN01370 

3CND1390 

SCN01400 


3CN01420 

iCNDUSO 

6CN01440 

3CN01450 


- 50 - 


SD 70-266 



space Division 

North American Rockwell 


<S3G,-630,63lt 

,-631,632, 

632, -632,633, 

633, -633,634, 

634, -634,63';, 

63>;,-635,636, 

636, -629,636, 

REM 

RFM fluid flow conductors (SYSTEM 3) 

REM 

637, -637,630, 

630,-630,639, 

639, -639,640, 

640, -640,641. 

641, -641,642, 

642, -642,643, 

643, -643,644, 

644, -644, 645, 

645, -645,646, 

646, -646,647, 

647, -647,640, 

640,-640,649, 

649, -649,650, 

650, -650,651, 

651, -651,652, 

652, -652,653, 

653, -653,654, 

654, -647,654, 

REM 

REM RADIATION CONDUCTOR - NUMBER 901 

-901,701,901,703,901,722,901,724,901 

725.902.727.902.746.902.748.902 

749.903. 751.903.770.903.772.903 

773.904.775.904.794.904.796.904 

797.905.799.905.018.905.620.905 

021.906.823.906.342.906.844.906 

045.907. 847.907.866.907.868.907 

069, 908, 871, 908, 890, 908, 892, 908,. 22982 7E-9 

REM 

REM RADIATION CONDUCTOR - NUMBER 902 
REM 

-902,702,901,723,901,726,902,747,902 

750,903,771,903,774,904,795,904 
798,905, 819,905,022.906,843,906 
846,9O7,867,9O7,07O,9O0,891,9Oe,.196995E-9 

REM 

REM RADIATION CONDUCTOR - NUMBER 903 
RE H 

-903,704,901,712,901,713,901,721,901 

728.902.736.902.737.902.745.902 

752.903.760.903.761.903.769.903 

776.904. 784.904.705.904.793.904 

800.905.800.905.809.905.017.905 

824.906.032.906.833.906.841.906 

840.907.856.907.857.907.865.907 

072, 908, 880, 908, 88 1,908,089, 908*. 523268E-9 

REM 

REM RADIAtlDN CONDUCTOR - NUMBER 904 
REM 

-904,705, ’ I, TCt , : I , ■ i 1 , ' - - , ■ I ,'/l i , ? 

729,902,732,9D2,7 3 5,9r'2,738,9" ' ,741,902,744,902 

753.903. 756.903.759.903.762.903.765.903.768.903 

777.904. 780.904. 783.904.786.904.789.904.792.904 

601.905.804.905.807.905.810.905.613.905.816.905 


tCN0l470 

SCND140O 

ACND149C 

SCN01500 


SCND152D 
$CND1530 
SCND154C 
tPANEL IACND1550 
tPANEL 2»CN01560 
IPANEL 3*CND1570 
iPANEL 4ACND158C 
iPANEL 54CN01590 
iPANEL 6«CND1600 
iPANEL 7tCN016l0 
iPANEL 8iCND1620 
tCNOl630 
iCND1640 
iCND1650 
iPANEL l-2iCND1660 
iPANEL 3-4iCNDl670 
iPANEL 5-&iCND168C 
iPANEL 7-8iCND1690 
iCN0l700 
iCNOlTlO 
iCNDl72C 
tPANEL ltCNOl730 
iPANEL 2iCND1740 
tPANEL 3iCNDl750 
iPANEL 4iCN01760 
tPANEL 5tCNDl770 
iPANEL 6iCNOl780 
iPANEL 7tCND1790 
iPANEL 8iCNDl800 
iCNOlOlO 
SCN01820 
:CN;U830 

it^’ikEL 1* ^ . 

tPANEL 2tCNO105O 
iPANEL 3tCN01660 
iPANEL 4iCND1870 
tPANEL 5iCNO109O 
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02 5,'>O6t 8a8,906,83l.'?06»834,<i06t837,906»S40t906 
099,007, 3 52, 907, 855 » 90 7, 85 8» 907, 661,90 7, 864» 907 
073,908,875,908,879,908,882,908,885,908 
888,908, .27C036E-9 

HEM 

PFM RAOUTION CONnUCTOR - NUM8ER 905 
R6'« 

-905,706,901,707,901,709,901,710,901 

715.901.715.901.718.901.719.901 

730.902.731.902.733.902.739.902 

739.902.740.902.78.2.902.763.902 

754.903.755.903.757.903.758.903 

743.903.764.903.766.903.757.903 


*P6NEL 6*CN0189C 
SPANEl 7tCN01900 
*CN01910 
tCND1920 
4CN01930 
SCN0194C 
*.CN0195C 
5PANEL 1SCN01960 
tPANEL UCND1970 
*PANEL 2KNO190O 
6PANEL 2ACND1990 
APANEL 3tCND20C0 
fPANEL 3<CND201C 


787.904.788.904 

802.905.803.905 

011.905.812.905 

826.906.827.906 

835.906.836.906 

850.907.851.907 

859.907.850.907 

874.908.875.908 

883.908.884.908 


,781,904, 

,790,904, 

,805,905, 

,814,905, 

,829,906, 

,838,906, 

,853,907, 

,862,907, 

,877,908, 

,886,908, 


782, 

791, 

006 , 

815, 

830, 

839, 

954, 

863, 

878, 

887, 


904 

904 

905 

905 

906 

906 

907 

907 

908 
908, 



PUT, 0.05 
4,10,5,27 »K5 IS 
, 0.0001 


tPANEL 43CN02020 
$PANEL 4$CN02030 
iPANEL 5$CN02040 
4PANEL 5SCN02050 
3PANEL 63CND2060 
5PANEL 63CN02070 
*PANEL 7$CNn2O0O 
SPANEL 7iCM02090 
iPANFL 6*CNn2100 
(PANEL 3»CN02110 
tCN0212C 
5CN02130 
(CNSOOlO 


NO PRESSURE CONO PEP SYSTEM (CNS0030 
$ WATER Gtvr.DL FLOW RATE 
( RVPASS OPTION 


359122E-9 


REM CONSTANTS FOR INCIDENT HEAT — PANELS 1-8 

12.. 13.. 14.. 15. 39.. 16. 83., 17, 82., 18, 40., 19, 

20,24 $ NUMBER OF NODES PER PANEL 

21,*374223E-3 t VOLUME OF FLUID CAPACITANCE NODES 

22,16 % NO. OF NODES PER SYSTEM WITH CAPACITANCE VALUES 

23,3.60 

REM -K22 IS NO. OF CONVECTION CONDUCTORS PER SYSTEM 

24,19 S NO. OF THERMAL FLUID FLOW CONDUCTORS PFR SYSTEM 

25.. 6944444,-2 * CONVERSION FROM PSF TO PSI 
REM CONSTANTS USED IN OUTPUT CALLS 

26. . 27. 16. 28. .29, STEMP,3C , 192, 31 ,701 , 32, . 33,N0 USFI 
REM MULTIPLYING FACTORS FOR PRESSURE CONDUCTORS 

501.255.91.502.108.265.503.108.265.504.30.57.505.108.265 
506,108.265,507,53.03,508,108.265,509,108.265,510,65.726 

511.108.265.512.108.265.513.49.71.514.108.265.515.108.265 

516.33.01.517.108.265.518.108.265.519.59.7692.520.108.265 
521,108.265,522,58.0113,523,108.265,524,108.265,525,36.0373 

526.55.2666.527. 148.1534.528.285.9462.529.108.265 
530,108.265,531,30.7636,532,108.265,533,108.265,534,67.6798 

535.108.265.536.108.265.537.49.71.538.108.265.539.108.265 

543.49.4174.541.108.265.542.108.265.543.66.5079.544.108.265 


tCN$ri060 

(CNS0070 

ICNS0080 

(CNS0090 

(CNSOlOO 

(CNSOtlO 

$CNS0120 

tCNSOlBO 

iCNSOUO 

SCNS0150 

SCNS0170 

$CNS0180 

(CNS0190 

*CNS0200 

$CNS0210 

tCNS0220 

$CNS0230 

(CNS0240 

iCNS0250 

(CNS0260 


545,108.265,546,43,9479,547,108.265,548,108.265,549,59.3824 SCNS0270 
550,108.265,551,108.265,552,47.0735,553,68.8417,554,159.4818 ICNS028C 
555,295.3176,556,108.265,557,108.265,558,46.7802,559,108.265 $CNS0290 


560,108.265,561,52.0539,562,108.265,563,108.265,564,50.296 $CNS0300 

565.106.265.566.108.265.567.65.2422.568.108.265.569.108.265 tCNS0310 

570.43.5577.571.106.265.572.108.265.573.44.4363.574.108.265 (CNS0320 

575, 108. 265, 576, 73. 8325, 577, 108. 265, 578, 108. 265, 579, 40. 9204 (CNS0330 

580,68.8717,581,158.5055 (CNS0340 

REM COMPRESSED PRESSURE CONDUCTORS (CNS0350 

582.. 583., 584,, 585,, 586,, 587,, 588,, 589,, 590,, 591,, 592,, 593, SCNS0360 

594.0. .595.7 5CNS0370 

598,150.0 (INITIAL SOLID NODE TEMP 


- 52 - 


SD 70-266 





Space Division 
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!i 99»192 t NH. nr SHLIO TEMPEt^ATUftE NnOFS tCNS 0390 

6 CC. 6.2 

REM CONSTANTS 601 TO 681 « 6 /WP AND ARE USEO IN RE NO. CONPUTATN SCNSO 610 
60 I, 66 .A 27 , 602 , 37 . 619 . 6 D 3 , 37 . 619 , 606 , 5 C. 097 . 6 C 5 . 37.619 tCNS 062 C 

606 . 37.61 9 . 607 , 50 . 087 , 608 , 37 . 61 9 , 609 , 37 . 619 , 610 . 50.087 trNS 06 30 

611 . 37 . 619 . 612 . 37 . 619 . 613 . 50 . 087 . 616 . 37 . 619 . 615 . 37.619 tCNSC 6<,0 

616.50.087.617.37.619.618.37.619.619.50.087.620.37.619 3CNS055C 

621.37.619.622.50.087.623.37.619.624.37.619.625.50.087 ^CNS0460 

626.53.087.627.50.087.623.66.427.629.37.619.630.37.619 $CNSC470 

631 . 50 . 097 . 632 . 37 . 619 . 633 . 37 . 619 . 634 . 50 . 087 . 635 . 37.619 $CNS 0480 

636 , 37 . 619 , 637 , 50 . 087 , 638 , 37 . 619 , 639 , 37 . 619 , 640 , 50.097 tCNS 0490 

641 . 37 . 619 . 642 . 37 . 619 . 643 . 50 . 087 . 644 . 37 . 619 . 645 . 37.619 *CNS 0500 

646 . 50 . 087 . 647 . 37 . 619 . 643 . 37 . 619 . 649 . 50 . 087 . 650 . 37.619 tCNSO 5 l 0 

651 , 37 , 619 , 652 , 50 . 037 , 653 , 50 . 037 , 654 , 50 . 067 , 655 , 66.427 tCNS 0520 

656 . 37 . 619 . 657 . 37 . 619 . 658 . 50 . 087 . 659 . 37 . 619 . 660 . 37.619 $CNSC 530 

661 . 50 . 037 . 662 . 37 . 61 9 . 663 . 3 7 . 61 9 . 664 . 50 . 087 . 665 . 37.619 ICNS 054 P 

666 , 37 , 619 , 667 , .08 - 68, 3 V. 6 r , -. 69 , j - y ,6 i'G, ' tC .SC'-Sl 

671.37.619.672.37.619.673.50.087.674.37.619.675.37.619 3CNS0560 

676 . 50 . 087 . 677 . 37 . 619 . 673 . 37 . 619 . 679 . 50 . 087 . 680 . 50.087 tCNS 0570 

681 . 50 . 037 . 682 . . 683 . . 634 , iCNS 0580 

701 . . 03031 1 . 702 . 0 . . 703 . 0 . . 704 . . 042633 . 705 . 0 . . 706 . 0 . *CNS 0590 

707 . . 042633 . 703 . 0 .. 709 . 0 .. 7 10 . . 042633 . 711 . 0 . . 712 . 0 . SCNS 0600 

713 . . 342633 . 714 . 0 .. 715 . 0 .. 716 . . 042633 . 717 . 0 . . 718 . 0 . tCNS 05 lO 

719 . . 342633 . 720 . 0 .. 72 1 . 0 .. 722 .. 042633 . 72 3 . 0 .. 724 . 0 , 5 CNS 0620 

725 . . 023644 . 726 . . 0023271 . 727 . . 04728 76 . 72 8 . . 080 3 11 3 CNS 0630 

729 . 0 . . 730 . 0 .. 73 1.. 0 42633 , 732 , 0 . , 733 , 0 ., 7 34 ,, 04263 3 , 735 , 0 . $CNS 064 C 

736.0. .737. .042633. 738.0. . 739.0. . 740 . . 0426 33 . 741 . 0. . 742 . 0. SCNS0650 

743 .. 042633 . 744 . 0 .. 745 . 0 .. 746 .. 042633 . 747 . 0 .. 748 . 0 . tCNS 0660 

749. . 042633. 750.0. . 75 1 . 0. . 752 . .02 3644. 75 3 . . 0C232 71 tCNS0670 

754. . 0472876.755. . 0803 1 1 . 756. 0. . 757. 0. . 758 04263 3 . 759. 0 , SCNS0680 

760 . 0 . . 761 . . 042633 . 762 . 0 . . 763 . 0 . . 764 . . 0426 33 . 765 ,C, , 766 , 0 . iCNS 0690 

767 .. 042633 . 768 . 0 .. 769 . 0 .. 770 .. 042633 . 77 1 . 0 .. 772 . 0 . tCNS 070 C 

77 3 .. 04263 3 . 774 . 0 .. 775 . 0 .. 776 .. 04263 ', 777 , 0 .. 77 P,C, $CNS 0710 

779 . . 023644 . 780 .. 002 32 71 . 78 1.. 0472876 $CNS 0720 

801. . 802. .803, t ELUIO HEAT LOSS PER SYSTEM $CNS0730 

804.. 0265823*(MP/41 FOR PANEL TUBES SCNS0740 

REM CONSTANTS FOR GR NUMBER « tWP/4)«(X/0l ACNS0750 

805.. 984934, 806, 4. 725, SOT,, 984934, 808, 4. 725 ICNS0T60 

809.. 984934. 8 10. 4. 72 5. 81 1.. 984934. 81 2. 4. 725 tCNSOTTO 

8 1 3.. 984934. 8 14. 4. 72 5. 81 5.. 984934. 8 16. 4. 725 $CNS0780 

817. . 984934.818.4.725.819. .984934.820.4.725 1CNS0790 

821. . 822. .823, iINVERSE FLOW RATE iCNSO0OO 

824. . 825. .826, S SYSTEM PRESSURE DROP SCNS0810 


827. 150.. 828. 150. .829. 150. 

830, 22.. 831, 22., 832, 22. 

833.. 73. 834.. 73. 835. .73 

836, 160., 837, 160., 838, 160. 

839, 85.8 

840, 767. 

841. 180.. 842. 180. .843. 180. 

844. 340.. 845. 340. .846.340. 

847. 420.. 848. 420. .849.420. 

850, -.60, 851,-. 60, 852, -.60 
853, 0.023 

854. 30. . 855.0. .856.0. 

857, O.Ol 

858. . 859. .660. 

861, 0. 

862, 30. 

863, 0. 

864,2,43E-3 
865,1.545E-2 
866,3. 76E-3 


S PREVIOUS TSRCE 
( MASS OF KOH 
% PERCENT KOH IN STACK 
S TCE 

% GAS CONSTANT - WATER 
i GAS CONSTANT H2 

t volumetric flow RATE 

$ STACK INLET TEMP 
t STACK OUTLET TEMPS 
t PC H20 PROD. IN STACK 
t H20 PRODUCTION RATE 
i CURRENT - AMPS 
» ACCEPTABLE ERROR FOR RATE 
* CELL VOLTAGE 
» NOT USED 
S MCP OF STACK 
S NOT USED 


*VOl ,FT**3,RA0 OUT TO REG COLO INLET 
iV0L,FT**3,REG COLO SIDE 
1VOL,FT**3,CONO GLYCOL SIDE 
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868,;!.^3lE-2 

867,2.26E-3 


tVJL,FT**3tC0ND OUTLET TO PEG. INLET 
3 ,FT*»3, SECONDARY REG HOT SIDE 

tVOL .FT**3.R£G. TO RAD INLET 

870. 160. . 871. 16C. ,872, 160. tTCE, PREVIOUS TIME STEP 

873, 1, 87A, 1 ,875, 1 MNTEGER POINTER FOR SFC RP VLV 

876. . 877. . 878, $ TEMP COND INLET 

878. . 5. 880. .'>,891, .5 $SECONOARY BYPASSED FRACTION 
882, 130. ,883, 130. ,888-, 13C. tSEC REG COLD OUTLET 

885. . ‘55 tBYPASS FRACTION CUTOFF VALUE 

886. . 887. .888, SINVERSE FLOM RATES 

FACTOR 

tREG COLD SIDE INLET 
$COND CLY OUT TEMP 
SSEC REG HOT SIDF INLET 
SSEC REG HOT SIDE OUTLET 

i INTEGERS USED FOR 
t INTEGERS 

908, 150.07,90'», 150. 07,910, 150.07 t COND CLY 


889,0.25 


*REG DELAY AOJ 



LOOet 

891, 

100- 

,892, 

100 


175*t 

894, 

175. 

,895, 

175 


llS.f 

897, 

175. 

,098, 

175 


lDO*f 

900, 

IOC. 

,901, 

100 

90?t 

t . .903,1 . 

,904 

t 1 9 


905 t 

,906, 

,907 

f 




DO LOOPS 
USED FOR BY PASS 
INLET TEMP 


911, 1. 

912,200 

913,199 

91A,20 

915,60 

916,30 


tCOND HEAT BAL TOLERANCE 
SARRAY IC 
SARRAY tC 
(ARRAY IC 
(ARRAY IC 
(ARRAY IC 


917, ,918, ,919, 

% 

POWER - WATTS 

920,77.29,921,77.29,922,77.29 

$ 

H20 SPECIFIC VOLUMF 

923, ,924, ,925, 

i 

PERCENT 8Y-PASS PRIMARY 

926, ,927, ,923, 

% 

STACK EXIT PRESSURE 

929, ,930, ,931, 

% 

PERCENT KOH 

932, ,933, .934, 

% 

H20 PRESSURE 

935, ,936, ,937, 

% 

ENTHALPY OF STEAM 

933, ,939, ,940, 

% 

INITIAL PRESSURE 

941, ,942, ,943, 

s 

INITIAL SPECIFIC VOLUME 

944, 70., 945, 70., 9- ,70, 

% 

AMBIENT TEMP OEG F 

947,0.00257 

i 

H2 PRODUCTION RATE 


948,5772.79 (LOWER HEATING VALUE, H2 

949,6824.42 (HIGHER HEATING VALUE, H2 

950,0.2 (ERROR FOR STACK INLFT 

(INITIAL CONO TEMP INCREMENT 
(INTEGER INDICATOR FOR COND IM8AL 

(ERROR FOR HEAT BALANCE 
(ERROR FOR STACK INLET TEMP 
(WATER CONDENSEO-LB/HR 
(STACK WATER STORAGE RATE 
( SEC BP VLV STAT POS 
( 

( 

( 

( 

( 

( 

( 


951, 5. 

952, ♦! 

953,1,0 

954.0. 2 

955, ,956, ,957, 

958. . 959. .960, 

961.0. 5.962.0.5.963.0.5 

964.530. 

965.530. 

966.0. 48 

967.0. 06 

960, IDO., 969, 100. ,970,100. 

971. 100.. 972. 100.. 973. 100, 

974.. 975.. 976, 


02 INLET TEMP - 
H2 INLET TEMP - 
02 FLOW RATE - 
H2 FLOW RATE - 
02 REACTANT TEMP 
H2 reactant TEMP TO STACK 
FLOW RATE FIX-UP 


OEGR 
OEGR 
LB/HR 
LB/HR 
TO STACK 


3ARRAY DATA 

1 ( DENSITY OF GLYCOL -WATER LBM/FT«*3 

-460,, 75, 09,-50,, 69, 54, -40,, 69, 4, -30,, 69, 25, -20,, 69.1 
- 1 0., 68. 93, O.t 68. 76, 10., 6ft. 58, 20., 68. 4, 30., 68. 2, 40., 68.1 

50. . 67.03. 60 • ,67, 63,70. ,67* 43 ,80*, 67*22, 90. ,67*, 100*, 66. 78 
UD.,66.56,12C.,66.34,130.,66.1,140.,65.85,150.,65.61 

163.. 65.35, t70.,65.t,180.,64.84,190.,64.57,200.,64.23 

250.. 62*79, 300., 61. 23, 350., 59. 61, END 

? ( CONDUCTIVITY OF GLYCOL-WATER BTU/HR*FTPF 

-460.,, 2 279, -50.,. 2279, -40.,. 2279, -30., ,2278, -20.,. 2277 


(CNS1050 

(ARROOlO 

(ARR0020 

1ARR0030 

SARR0040 

(ARR0050 

(ARR0C60 

(ARROOTO 

(ARR0080 

(ARROC90 

(ARROlOO 
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-n,,,2?75,C.,.2a73aO.,.22 70,20.,.226f>,30.,.2?6?,^0.».225fr 
5C. , .2252,60. , .2246,70, , , 2 240 t 30 . , , 2233 ,90. ,.2225, IOC. ,.2217 

1 10.. . 2208. 1 20,,. 2 199, no.,. 2 189, 140.,. 2 178, 150.,. 2 1 67 

160.. . 2156. 1 70... 21 43. 180. ..21 30. 190.. .211 7. 200.. .2103 

340.. . 1844, FND 

3 % SPeCIFIC HEAT OF GLYCOl -MATER 3TU/t,3H*F 
-460. , .033,-50. , .606,-40. , .618,-30. , ,631,-20, , .643 
-10. ,.655,0. , .666,10. , .677,20, , , 688 , 30, , . 693 ,40, , . 708 
5C.,. 718,60., .728,70,, . 737 , 80 . , . 746 , 90. , , 754, 1 00 . , , 762 

1 13. . . 770. 120. . . Y77, 130., .784, 140. , . 70l , 1 50. , . 797, 1 60 . , . 803 

173. . .809. 180. . .814. 190. . .819.200.. .824.230. 84 
270. , .36,310. , .878, HMD 

4 i VISCOSITY OF GLYCOL-WATER E-3 LRH/FT*SEC 

-460.,50DOOO.,-50.,245.,-40.,135,,-30.,80.,-20.,51.9,-10. 

34.. 0.. 24. 5. 10.. 16. 5. 20. . 12. 2. 30.. 9. 3. 40.. 7. 3. 5C., 5. 75 

63. . 4.65. 70., 3. 75,80., 3.05,90, ,2.58, 100., 2.08, 1 10. ,! .87 

120. . 1 .61 1 130. . 1.4. 140.. 1. 22. 1 50. ,1 . 11, 160. ,.96, 170. , .86 

180.. . 77. 1 90... 69. 200... 625. 2 50... 406. 300... 303 

350.. .216, END 
5,SPACE,28,END 
6,SPACF,28,END 
7, SPACE, 28, END 
8, SPACE, 81, END 
9, SPACE, 54, END 

10 tBYPASSED FRACTION — SECONDARY SP VLV- ASCEND 1 NG TCE 

157.. 0.0, 157. 25, 0.1 65, 157. 5, 0. 245, I 57, 75, 0 , 3 , 1 58. , 0. 34 

1 58. 25. 0. 3 0.1 58. 5. 0.41. 158. 75. 0.43. 1 59.. 0.445. 1 59. 2 5. 0.46 
I 59.5,0.47, 1 59. 75, 0.474, 160,, 0,48, 160.2 5,0.485,160,5,0.49 

160.75.0. 493, 16 1 . , 0 .498 , 16 1 , 2 5, 0 . 5 , 1 6 1 . 5 , 0 . 5 1 , 161 . 75 , 0. 52 

1 62.. 0. 534. 162. 25. 0.555. 162. 5. 0.58. 162. 75. 0.61 5.1 63., C. 65 

163.25.0. 71. 163.5.0.8.163.75.0.935. 164. .1 .0, END 

11 ^BYPASSED FRACTION — SECONDARY BP VL V-OESCENO TNG TCE 
I 54,, 0,0, 154.25,0.12,154.5,0. 195, 154.75,0.255, 155. ,0.31 
155.2 5,0.365, 1 55 . 5 , 0. 4, 1 55. 75 ,0. 425 , 1 56. , 0.44 , 1 56 . 2 5,0. 457 

156.5.0. 467.156.75.0.475.157.0.0.48. 157.25.0.485.157.5.0.488 

157. 75. 0. 492. 1 58.. 0.5. 158. 25. 0.506. 1 58. 5. 0.5 1 2. 158. 7 5, 0.522 

159.. 0.53 5. 159. 2 5. 0.554. 1 59, 5, 0.58, 159. 75, 0,62, 160. ,0,665 

160.25.3.74.160.5.0. 88. 160. 75.0. 94. 16 1 . *1 .0 , END 

12.. 145833,, 125000, .145833,, 33203 1, ,171 347,. 227875,. 227075 
.171347,. 2270 75,. 2278 75,. 1 71347,. 332031,, 332031,. 171347 
.2278 75,. 227875,, 171 347,. 2278 75.. 227875,. 171 347,. 332031 
. 145833,. 125000, .1458 33, END 

REM THE FOLLOWING ARRAYS PROVIDE FDR VARIABLE INCIDENT HEAT FLUX 

13 » INCIDENT HEAT FOR PANEL I - BTU/HR FT**2 
0.00,07.09,0.05, 15.67,0.10,23.90,0.15,31.45,0.20,38.00 
0.25,43,30,0.30,47.13,0.35,49.33,0,40,49.81,0.45,48.55 
0.50,45.60,0.55,41.09,0.60,35.20,0.65,28.16,0.70,20.27 
0.75, 11.84,0.80,07.09, 1.55, 07, 09, END 

14 $ INCIDENT HEAT FDR PANEL 2 - BTU/HR FT**2 
0.00,00.27,0.05,15.53,0.10, 30.26,0.15,43.73,0.20,55.43 
0.25,64. 88,0.30,71.71,0.3 5,75.63,0.40,76.48,0.45,74.24 
0.50,68,98,0,55,60.94,0.60,50.42,0,65,37.06,0.70,23.77 
0.75,08.74,0.80,00.27, 1.55, 00. 27, END 

15 % INCIDENT HEAT FOR PANEL 3 - BTU/HR FT**2 
0.00,00.00,0.05,16.82*0.10,32.96,0.15,47.75,0.20,60.61 
0.25,71.00,0.30,78.50,0.35,32.81,0.40,83.75,0.45,61.28 
0.50,75.51,0.55,66.66,0.60,55. 11,0.65,41.31,0.70,25.83 
0.75,09.31,0.80,00.00* 1 . 55, 00 *00 , END 

16 i INCIDENT HEAT FOR PANEL 4 - BTU/HR FT**2 
0.00,06.91,0.05,16.24,0.10,23.28,0.15,20.74,0.20,35.34 
0.25,39.88,0.30,43.15,0.35,45,03,0.40,45.44,0.45,44.36 
0.50,41.04,0.55, 37, 99, 0.60, 32, 94, 0.65, 26. 93, 0.70,20. 18 
0.75, 12. 97, 0.80, 08, 9 1, 1.55, 00. 91, END 
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tARROl 10 
lARftOl/C 
TARROnO 
SARR014C 
tARRClSC 
t ARR016C 
$ARRC170 
$ARR0t30 
tAWROno 
4ARR0200 
IARR0210 
tARR0211 
tARRC220 
tARR0230 
tARRC24C 
tARRC250 
SARR0260 
$ARR0?70 
tARR027l 
$ARR02RO 
SARR0290 
$APR0300 
SARR0310 
SARR0320 


$ARR035C 

♦ARR0360 

SARR0370 

tARR0380 

$ARR0390 

IARR0400 

tPANELl 

SPANELl 

iPANFLl 

iPANELl 

tARR0420 

5PANEL2 

»PANEL2 

IPANEL2 

iPANEL2 

8ARR0440 

8PANEL3 

SPANEL3 

$PANEL3 

I PANEL 3 

5ARR0460 

SPANEL4 

5PANEL4 

(PANEL4 

$PANEL4 
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1? % INCIDENT HEAT FOR PANEL 5 - STU/HR FT**? 

0.3C.32.37,O.P5,35.AUO.l0.39.3?»0.15t4l.OOtO.?OtA3.3a 
C. 25t 45. 2C, 0. 30, 46. 56, C. 3 5, 47. 34, 0.40*4 7. 5 1, C. 45, AT, 06 
0.50, 46. 02, 0.55, 44.41,0.60,42,33,0.65,39.8 3,0.70,37.04 
0.75, 34.05,0.80,35.82,0. 85,43.35,0.90,32.37, i .45,32.37 
1.5C, 40,76, 1.55, 32. 37, END 

18 i INCIDENT MEAT FOP PANEL 6 - BTU/HR FT**2 

3.00,49. 71,0.05,54.38,0. 10, 58.86,0. 15, 62.67, C. 20, 66. 54 
0.25, 69. 42, 0,30, 71. 50, 0.35,72.70,0,41', 72. 96, 0.45, 72, 2 8 
0.53,70.67,0.55,68.22,0.60,65.01,0.65,61.18,0.70,56.88 
0.75,52.29,0.80,56,54,0.85, 71,47,0.90,49.71,1.45,49.7] 
1 .53,66.34, 1. 55, 49.7 1, ENO 

19 S INCIDENT HEAT FDR PANEL 7 - BTU/HR FT**2 


0.30,56. 74, 0.05, 6 2.08, 0.1 0,67, 1 9, 0.1 5, 7 1.6 8, 0.20,75. 95 
0.2 5, 79, 2 5, 0.30, 8 1.62, 0.35, 82. 99, 0.40, 8 3. 29, 0.45, 32. 50 
0,50,30.67,0.55,77,87,0.60,74,21,0.65,69.84,0.70,64,93 
0.75,59.69,0.80,64.60,0.85,81.75,0.90,56.74,1.45,56.74 
1.50, 75. 86, 1.55, 56. 74, END 

20 » INC IDENT. HEAT FOR PANEL 8 - BTU/HR FT**2 

0.00,31.81,0.05,34.80,0.10,37.66,0.15,40.29,0.20,42.57 
0.25,44.42,0.30,45.75,0.35,46.52,0.40,46.68,0.45,46.25 
0.50,45.22,0.55,43.65,0.60,41.60,0.65.39.15,0.70,36.40 
0.75,33.46,0.80,35. 13, 0 . 85, 42 . 39 ,0 .90, 3 1 ,8 I , 1 . 45 , 3 1 . 8 1 
1.50,39.90, 1.55, 3 1.8 1, END 


21, SPACE, 56, END 

22, SPACE, 81, END * REYNOLDS NUMBER 
23, SPACE, 81, END 

24 SRELATlve ENTHALPY OF GLYCOL-HATER, BTU/LflH 
-460., 0.0, -50., 131. 0,-40., 137.12,-30., 143. 36, -20., 149. 73 
-10. , 156.22,0.0. 162.83, 10. , 169.54,20. . 1 76.37,30. , 183. 3 
4C., 190. 33, 50.. 197. 46, 60. , 204. 69 , TO . , 21 2. 0 1 * 80. , 2 1 9. 43 

90.. 226. 93. 10C., 2 34. 5 1, 110., 242. 17, 120., 249. 9, 130., 25 7. 71 

140.. 265. 58. 150.. 273. 52. 160.. 281. 52. 170.. 289. 58. 180.. 297. 7 
190, *305. 86, 200., 31 4. 08, 2 30., 339. 04, 2 70., 373. 04, 3 to., 407. 8 
END 

25, SPACE, 48, END 

26 A on LOOP INDICES - INITIAL VALUE 

0. , l.,6.2, l.tEND 

27 * DO LOOP INDICES - TEST VALUE 

0., 1* ,6.2, l.,END 

28 $ 00 LOOP INDICES - INCREMENT 

0. ,1 *,6.2t 1. tENO 

29, SPACE. 10, END I NOT USED 

30, SPACE, 10, END % NOT USED 

31, SPACE, 10, END $ NOT USED 

32, SPACE, 10, END A NOT USED 

33, SPACE, 10, END A NOT USED 

34, SPACE, 10, END A NOT USED 

35 A BYPASS OPTION - SYSTEM I 

0. ,0. ,3*1 ,0.,3* 101, 1*, 6. 2,1., END 

36 A BYPASS QPTI? j - SYSTEM 2 

0. ,0. ,6.2,0. ,END 

37 A BYPASS OPTr’^; - SYSTEM 3 

0 . , 0 . , 6 * 2 , 0 . , END 

38 A HEAI TRANSFER COEFFICIENT - AREA TERM 

11. 11 .50. 20.. 0. .20. .40. .60. .80.. 100. . 120..14C* ,160* ,180. 

120.. 209..243..265..289..3t8.,345.,3T7.,0.,0.,0.,0. 

130*, 205. ,238. ,256* ,281*, 308*, 333. ,362. ,404* ,0. ,0. ,0. 

140.. 202.. 234.. 253.. 275.. 300*, 323., 350., 387., 0.,0.,0. 

150*, 199., 230., 248., 269., 293., 31 5., 339., 374. ,408., 0.,0. 

160.. 197.. 227.. 244.. 264.. 287.. 307.. 330.. 363.. 395.. 0..0. 

170.. 195.. 224*, 241., 260., 782., 301., 323., 353., 383., 41 7., 0. 

180.. 193.. 221.. 237.. 256.. 276.. 295.. 315. *344., 371. ,404.,0. 
t90.,190.,2l8.,234.,25l.,272.,289.,306.,335.,36l.,39t.,4lT 


APANEL5 

APANEL5 

APANEL? 

*PANEL5 

tPANFL5 

AARROSOO 

APANEL6 

APANbLS 

APANFL6 

APANE16 

APANEL6 

A ARR0520 

AP4NFL7 

APftNEL? 

APANEL7 

APANEL7 

APANEL7 

AARR0540 

APANFL8 

APANEL8 

APANEL8 

APANFL8 

SPANEL8 

AARR0560 

AARR057P 

AARR0580 


AARR06C0 


ATVAOOOl 

$rvA 0002 

ATVA0003 

ATVA0004 

ATVA0005 

ATVA0006 

ATVA0007 

ATVAOOOS 

ATVA0009 
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2D3, * t 89. ,216. ,249., 268. 

213. . 189. .216. .231. .248. .268. 

223.. 1 . 87., 2 14., 2 29., 245., 264. 
11 , 1 1 ,55. ,'*20.,0. ,20. ,40. ,60. 

123.. 223.. 263.. 287.. 3 12.. 342. 

133.. 223.. 258. .280.. 304.. 332. 

143.. 2 16.. 254.. 275.. 297.. 32 3, 
153. ,213. ,250., 270., 291., 316. 

163.. 2 11.. 246.. 265.. 285.. 309. 

170.. 208.. 243.. 262.. 281.. 304. 

18 3.. 206.. 240.. 258.. 2 76.. 298. 

193.. 203.. 2 37.. 254.. 2 72.. 293. 

20 3.. 20 2.. 2 34.. 251. .269.. 289. 

213.. 202.. 234.. 251.. 269.. 289. 

223.. 203.. 232.. 248.. 265.. 285. 
11 ,11 ,63. ,'”20.,0. ,20. ,40. ,60* 

123.. 240.. 2 85.. 310.. 334.. 365. 

133.. 236.. 279.. 303.. 32 5.. 354, 

143.. 2 33.. 274.. 296.. 31 8.. 345. 

153.. 230.. 2 70.. 291.. 311.. 3 37. 

163.. 227.. 266.. 286.. 336.. 330. 

173.. 2 25.. 263.. 282.. 301.. 324, 

18 3.. 222.. 2 59.. 2 78.. 296.. 319. 

193.. 223.. 256.. 2 74.. 292.. 3 1 3. 

23 3.. 21 8.. 253.. 271.. 288.. 309. 

2 13.. 2 18.. 253.. 271.. 288.. 309. 

223.. 2 16.. 252.. 268.. 285.. 305. 
11, It ,65. ,^20*,0* ,20., 40. ,60* 

123.. 256.. 305.. 331.. 354.. 385. 

1 33.. 252. .299.. 324.. 345.. 374. 

143.. 240.. 294.. 3 1.. ,338., 364. 

153.. 245.. 289.. 311.. 331.. 356. 

163.. 242.. 285.. 307.. 3 25.. 349. 

173.. 243.. 281.. 302.. 320.. 343. 

18 3.. 237.. 278.. 298.. 31 5.. 33 7. 

193.. 234.. 2 74.. 294.. 310.. 331. 

233.. 233.. 272.. 291.. 30 7.. 327. 

2 1 3.. 232.. 272.. 291.. 30 6.. 327. 

223.. 231.. 269.. 287.. 33 3.. 323. 

11. 11. 73.. -20., 0., 20. *40., 60. 

1 23.. 266.. 322.. 350.. 373.. 403. 

133.. 262.. 316.. 343.. 364.. 392. 

140 . . 258. . 310 . . 336. . 356. . 382. 

153.. 255.. 306.. 330.. 349.. 374. 

163.. 2 52.. 301.. 325.. 343.. 367. 

173.. 249.. 298.. 320.. 337.. 360. 

183.. 247., 294., 316,, 332., 354. 

193.. 244.. 290.. 311.. 327.. 348. 

203.. 242. *287., 308., 324., 344. 

2 13.. 242.. 287.. 308.. 323.. 344. 

223.. 240.. 284.. 305.. 320.. 339. 

11.11 .75.. -20.,0.,20.,40.,60. 

123.. 266.. 334.. 366.. 389.. 420. 

133.. 262.. 328.. 358.. 379.. 408. 

143.. 258.. 322.. 351.. 371.. 398. 

153.. 255.. 317.. 345.. 364.. 390. 

163.. 252.. 3 13.. 340.. 358.. 383. 

173. . 253.. 309.. 335.. 353.. 376. 

1 83.. 247.. 305.. 330.. 347.. 370. 

193 . . 245. . 301 . . 326. . 342. * 364. 

203.. 243.. 299.. 322.. 338.. 359. 

2 13.. 243.. 299.. 322.. 338.. 359. 

223.. 241.. 296.. 319.. 334.. 354. 


,285. ,303. 
,284. ,303. 
,200. ,290. 
,00. ,100*, 
,373. ,409. 
,360. ,393. 
, 349. ,380, 
, 340, ,369. 
,332. ,359. 
,325. ,351. 
,319. ,343, 
,313. ,335. 
,308. ,330. 
, 308. ,330. 
,303. ,324. 
, 80. , 100 . , 
,398. ,438. 
,384. ,421. 
,373. ,407. 
,363. ,395. 
.355. ,385. 
,348. ,376. 
,341. ,368. 
*335. ,360. 
,330. ,354. 
,330. ,354. 
,325. ,340. 
, 80 * ,100*, 
,421. *462. 
,407., 445. 
,395. ,430. 
,385. ,418. 
,377. ,408. 
,369. ,399, 
,362. ,390. 
,355. ,382. 
.350. ,376. 
,350. ,376. 
,345. ,369. 
, 60 * , 1 OC * , 
,441. ,484. 
,427. ,467. 
,415. ,452. 
,405., 439. 
,396. ,428, 
,389. ,419. 
,361. ,410. 
,374, ,401. 
,369. ,395. 
,366. ,395. 
,363. .389. 
, 60 * , 1 00 . , 
,460. ,504. 
,446, ,486. 
,433., 470, 
,423. ,458. 
,414., 447. 
,406. ,437. 
,399. ,428. 
,391. ,419. 
,386. ,413. 
,386. ,413. 
,380. ,406. 


, 329. ,354. ,381. 
, 329, , 354. ,381 . 
,323. ,346. ,372. 
12C., 140. ,160. , 

,0.t0.,0«,0. 

,434. ,0* ,0. ,0. 
,417. ,0* ,0. ,C. 
,403. ,438. ,0. ,0 
,391.,423.,C.,0 
, 30 I, ,411, ,449. 
,371. ,399. ,435. 
,362. ,383. ,421. 
,356. ,301. ,412. 
,355. ,380. ,411. 
,348. ,372. ,402, 
120. ,140. ,160. , 
,0* ,0*,0. ,0. 
,462*.0. ,0* ,0* 
,444. ,0*,0* ,0. 
,429. ,467. ,0. ,C 
,417.,452.,0.,C 
,407. ,439. ,478. 
,396. ,426. ,463, 
,387. ,415. ,449. 
,380. ,407, ,439, 
,380. ,406. ,439. 
,372. ,398. ,429. 

120. . 140. . 160. , 
,C* ,0. ,0*,0* 
,488, ,0. ,0* ,0. 
,470. ,0* ,0* ,0. 
,454.,494.,0.,0 
f 442 * , 4 77 * 1 0 . , 0 
,431. ,465, ,505. 
,420, ,452. ,489. 
,410. .440. ,474. 
,403. .432. ,464. 
,403. *431. ,464. 
,395. ,423, ,454. 

120. . 140. .160., 
,0*,0*,0*,0* 
,512. ,0* ,0. ,0* 

, 49 3*,0*,0. ,b* 
,478.,52l.,0.,0 
,465.,505.,0.,C 
,453. ,491, ,529, 
,442, ,477. ,513, 
,432. ,465. ,498. 
,425. ,456. ,488. 
,424. ,456. ,487. 
,417. ,447. ,477, 

120. . 140. .160., 
,0.,0*,0* ,0* 
*535* , 0. ,0* ,0* 
,515. ,0* ,0.,0V 
,499. ,546* ,C* ,0 
,486*,529.,0.,0 
,474. ,515. ,552. 
,463. ,501. ,536. 
,452. ,488. ,520. 
,445. ,479. ,510. 
,445. ,479. ,509. 
,437. ,469. ,498. 



$TV0OO1O 


tTVftOCU 

, 395 * 

tTVft00l2 

180* 

$TVA0013 

»TVftOOl4 

tTVA0O15 

«TVa00l6 


S7V4001 7 


ST VACO 10 

,0* 

1TVA0019 

,0* 

STVA0020 


$TV40C2t 


STVA0022 

,440, 

STVA0023 

,429* 

STVACC24 

UO* 

STVA0025 

STVA0026 

STVA0027 

STVAOO20 


STVA0029 


STVA0030 

,0. 

STVA0031 

,0* 

trVA0032 

,483* 

SrvA0033 

,472* 

ITVA0034 

,470. 

STVA0035 

,459. 

$TVAC036 

lao* 

STVA0C37 

STVA0038 

STVA0039 

STVA004C 

* 

STVA0C41 


STVA004? 

,0. 

STVA0043 

,0. 

STVA0044 

,511. 

STVA0045 

,499. 

STVA0046 

,498. 

STVA0047 

,486. 

$TV*0048 

180. 

STVA0049 

1TVA0050 

iTVA005l 

STVA0052 

* 

STVA0053 


STVA0054 

,0. 

STVA0055 

tO. 

STVA0056 

,536. 

STVA0057 

,524. 

STVA0058 

,523. 

STVA0059 

,510. 

STVA0C60 

180. 

STVA0061 

STVA0062 

STVA0063 

STVA0064 


STVA0065 


STVA0066 

.0. 

STVA006T 

,0. 

STVA0068 

,559. 

STVA0069 

,546. 

STVA0070 

,546. 

STVA0071 

,533. 

STVA0072 
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)U 1 1 ,2 0.,40,»60. flO. tlOO.t 120.» UO. , 16C. 

123, ,2 54. ,33T. ,'♦77. ,522. ,C.,C. .0. tO, 

1 3D. ,250 . ,33 3., 36fl. ,392.,42 3. ,463. ,504. , 555. ,0. ,0. ,0. 

14D. ,247. ,327. ,361. ,38 4., 4 13. ,450. , 488. , 536. , C , , 0 . , C . 

15D, ,244. , 322. ,355. ,377. ,404. ,439. ,475. , 519. ,568. ,0. ,C. 
16D., 242., 318., 350., 370., 397., 430. ,464. , 506. , 5 52 . , 0 . , 0 . 

X7D. ,239. ,314. , 345. , 365., 300. ,422. ,455. ,494. ,537, ,574. ,0. 
18 D. , 2 37, ,310,,34n. , 359 . , 384. , 41 5. , 445 . , 483 . , 5 2 3 , ,557. ,0, 

193. . 235. . 307. . 336. . 354. . 378. . 407. . 436. . 4 72. . 5 1 0 . .542. ,581. 

203. . 233. . 304. .333. . 350 . . 373 . . 402 . . 430. . 464. . 501 . . 53 1 . . 568 . 
21D.,233.,3C4.,332., 350 . , 373 . , 40 1 , , 430 , , 464. , 500 . , 53C, ,567. 
22D., 2 31 ., 30 U, 3 29., 346., 368., 396., 42 3., 456., 49 U, 5 1 9., 554. 
11, n ,85. ,-20., 0, ,20, ,40., 60. ,80. , 100. , 120. , 140. ,160. , 180. 
I2D. ,230. ,335* ,380* ,41U.,647. ,493* ,540* ,^« ,0* ,0* ,0* 

1 30.. 2 27.. 329.. 3 72.. 401.. 436.. 478.. 5 21.. 574.. 0..0..0. 

1 40.. 224. . 324. . 365. . 393 . . 426. . 465 . .506. .555..0..0..0. 

1 50.. 222.. 3 19.. 359.. 38 5.. 4 17.. 454.. 492.. 5 38.. 5 89.. 0..0. 

160., 220., 3 15., 354., 379., 4 10., 445., 4 0U, 524., 5 72., 0.,0. 

170. . 218. .312. .350. .374.. 403. .43 7. .471. .512. .557. .596. .0. 

180., 2 16,, 308., 345., 368,, 397., 429., 462,, 501., 5 43., 579., 0. 

190. . 215. . 305. . 341. . 363 . . 390. . 422 . . 453 . . 490 . . 529. . 563.. 603. 

200. . 213. .3C2. ,337, , 359., 386. , 41 6. , 447. , 48 2. , 520 . , 552 . , 590. 

2 10.. 2 1 3.. 302.. 3 37.. 359., 38 5., 4 16., 446., 48 2., 5 19., 552., 5 89. 

220.. 2 1 2., 299., 334., 355., 381. ,410., 439., 473., 5 1C., 540., 580. 
11, 1 1 ,90. ,-20,,C. ,2G.,40.,60. ,80. ,100., 120. ,140. , 160. ,180. 

120. . 199.. 322. .377. . 41 4. . 459. . 507. . 558. . 0. . 0. , G . , C . 

130.. 197.. 31 7.. 3 70.. 405.. 447.. 492.. 539.. 592.. 0..0..0. 

140.. 194. .312.. 363.. 397.. 437.. 479.. 523.. 572.. 0..0..0. 

150.. 193. . 308.. 358. .390.. 428. .46a.,5l0.,556.,606.,0.,0. 

160., 1 91., 30 5,, 35 3., 384., 421., 459., 499., 542., 589., 0.,0. 

170. . 190.. 301.. 348. .379.. 41 4. .451. .489. . 530. . 574. . 61 8 . .0 . 

180. . 189. .298.. 344. .374.. 408. .443. .479. . 5 1 8 . . 56C . , 601 . ,0 . 

190., 187, ,295. ,340, , 369. , 40 1 , , 435 . , 470. , 507. , 546 . ,585., 626. 

200.. 186..293. . 33T. ,365. , 397, , 43C . , 463. ,499. , 5 37. , 574 , , 6 1 3. 

210*, 186. ,293*, 36 *,397*,4_‘.. *,462^ * j*,5>6.r »* ! 

220.. 185.. 290.. 3 33.. 361.. 392.. 424.. 456.. 490.. 527.. 562.. 598. 
END 

39 $ SPECIFIC VOLUt^E OF STEAW - FT**3/L6H 


no 

.OOt 

157 

.34,132. 

00, 

149 

.66,134 

.00, 

142 

*42 

,136 

*00, 

135 

.58 

138 

• OOt 

129 

. 12,140. 

00, 

123 

.01.142 

.00, 

117 

*23 

,144 

.00, 

111 

.77 

U6 

.00* 

106 

.6C, 148. 

00, 

toi 

.71,150 

*00, 

97 

.07 

,152 

.00, 

92 

.60 

154 

*00, 

88 

.52* 156. 

00, 

04 

.58,158 

*00, 

00 

*84 

,160 

.00, 

77 

,29 

162 

*00, 

73 

.92, 164. 

00, 

70 

.73,166 

.00, 

67 

.69 

, 160 

.00, 

64 

.6C 

17C 

*00, 

62 

.06,172. 

00, 

59 

.45,174 

.00, 

56 

.97 

,176 

*00, 

54 

*61 

17$ 

*00, 

52 

.37, 100, 

CO, 

50 

.23. 182 

.00, 

40 

.19 

,184 

*00, 

46 

.25 

186 

.00, 

44 

,4C, 188. 

00, 

42 

.64,190 

*00, 

40 

*96 

,192 

.00, 

39 

*36 

194 

*00, 

37 

.83, 196. 

00, 

36 

.37,198 

.00, 

34 

.97 

,20C 

.00, 

33 

.64 

202 

*00, 

32 

.37,204. 

00, 

31 

.15.206 

*oc. 

29 

*99 

,208 

*00, 

20 

.88 

40 


$ 

SATURATl 

ON 

PR6 

SSUWE tlF STEAM 

- 

P$ 

!A 



13C 

.00, 

2 

.22, 132. 

00, 

2 

.34,134 

.00, 

2 

.47 

, 136 

*00, 

2 

.60 

138 

.00, 

2 

,74,140, 

00, 

2 

*89,142 

*00, 

3 

.04 

,144 

.00, 

3 

,20 

146 

*00, 

3 

.36, 148. 

00, 

3 

*54,150 

.00, 

3 

*72 

*152 

*00, 

3 

*91 

154 

*co. 

4 

.10, 156. 

00, 

4 

*31,158 

.CO, 

4 

*52 

,160 

*00, 

4 

*74 

162 

*00, 

4 

.97,164. 

00, 

5 

*21,166 

*00, 

5 

.46 

,168 

*00, 

5 

.72 

170 

*00, 

5 

.99, 172. 

00, 

6 

.72,174 

.00, 

6 

.56 

,176 

.00, 

6 

*87 

178 

*00, 

7 

.18,180. 

00, 

7 

*51, 102 

.00, 

7 

.05 

,164 

*00, 

0 

*20 

106 

.00, 

8 

.57, 188. 

00, 

0 

*95,190 

.00, 

9 

*34 

,192 

*oc. 

9 

*75 

194 

*00, 

10 

.17,196. 

00, 

10 

*61,190 

.00, 

11 

.06 

,200 

*00, 

11 

*53 

202 

* 00 , 

12 

.01,204. 

00, 

12 

.51,206 

*00, 

13 

.03 

,200 

*00, 

13 

.57 

^1 


% 

CURRENT 

- 

AMPS - 

SYSTEM 

1 






0. ,3C. ,3. l,30.,3. 1C 1, 25., 6. 2, 25. , END 
42 $ CURRENT - ANPS - StSTEH ? 

0. *0*, 6*2,0*, END 


trva0073 
$rVA0074 
tTVA0nT5 
tTVA0f!76 
tTVAC077 
tTVA007A 
tr VA0079 
3TVA008C' 
MV60C81 
3TVA0082 
tTva0D83 
*T VA0C84 
1.TVAC085 
STV^>0084 
tTvaC087 

trvaoo88 

tTVA0089 

tTva0090 

$TVA0C9l 

STVA0C92 

3TVA0093 

»TVA0094 

4TVA0095 

tTVA0096 

TTVA0097 

tTVA0098 

tTVAC099 

5TVA0100 

3TVA0101 

5TVA01C2 

tTvaoi03 

iTVACl04 
tTVAOlCS 
5TVA0106 
»Tv OlO: 
ATVA0108 
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4^ t CURRENT - ANPS - SYSTEM "i 

0. »C. t6.2*0. f END 


44 A PARTICAL 

PRESSURE OF 

WATER 

IN STACK - 


3. 0.55, 0.60 

, C * 65 

, 0.70 

* 0.75 

* 0 . 80 * 

0.05, 

C.90 

360.0, 27.0, 15.0, 

8.0, 

3*7t 

L.6t 

c.o* 

0*0, 

0*0 

33C.0, 36. C, 21.0, 

U.4, 

5*6, 

2.4* 

l.O, 

Q.O* 

c.o 

400. D, 47.0, 29.0, 

16.0, 

3*0, 

3.4* 

1*4, 

0.0* 

0.0 

420.0, 51.0, 38.0, 

21.4, 

10*8* 

■ 4 * 8 * 

1.9, 

O.C, 

0.0 

440,0, 79.0, 51.0, 

29.0, 

I4*6f 

6.6* 

2.7* 

l.C, 

c.o 

^6Dt0«lO2*Ot 66*0* 

38.0, 

19.0, 

8.9* 

3*6, 

1.3, 

0.0 


50.0, 

Z5.0* 

11.9* 

4.8* 

1.7* 

c.o 

5CC.0» 160<0» 104, Ct 

65.0, 

32*5, 

15.6* 

6*3, 

2.2, 

0.0 

52C*OtZOO.O,134,Ct 

83.0, 

42*0, 

20.6, 

0.4* 

3.0, 

1.2 

54C.Ct250*0, i70*0» 

104.0, 

55.0* 

27.5* 

11*3* 

4.2* 

1.6 

560,Cf3lO,0,2lO.Ot 

130.0, 

70.0, 

36.0* 

15.0, 

5.9, 

2*1 

533,0*380. 0* ?60-0, 

16C.0, 

69,0* 

46,0* 

20.0, 

0.0, 

2.3 

600*0*470*0, 32C^0* 

195. C, 

U3.0, 

60.0* 

27.0* 

10.7, 

3. 7, END 


45 t ENTHALPY OF SUPERHEATED STEAM - BTU/LBM 


5, 

126.08 

200.0* 

1150.0, 

220.0* 

1159.1, 

240.0, 

1168.2, 

260*0* 

1177. 3, 

260.0* 

1106.5* 

300*0* 

1195.6* 

320.0* 

1204.8* 

340.0* 

1213.9, 

36D.0* 

1223.1, 

380.0* 

1232.3* 

4DC.0* 

1241.6, 


4A, SPACE. l»END 
47, SPACE. 30. END 
4S. SPACE. 30. END 
4'!. SPACE. 60. END 
SO. SPACE. 199, END 
51. SPACE, 20, END 
52, SPACE. 20. END 
53, SPACE, 20, END 
54, SPACE, 200. END 
55. SPACE. 200. END 
56, SPACE, 200. END 
57. SPACE. 200, END 
5A, SPACE, 200, END 
59, SPACE, 200. END 
60, SPACE. 40, END 
61, SPACE, 40.EN0 
62, SPACE. 40. END 
63. SPACE, 200, END 
64. SPACE, 200, END 
65, SPACE. 200. END 


152. 97» 

170*06 

1149.2, 

1148*3* 

1158.5, 

1157.7* 

1167.6, 

1167.0, 

1176.8, 

1176.3* 

1186.0, 

1185.5* 

1195.2, 

1194.7* 

1204.4, 

1204*0* 

1213.6, 

1213*2* 

1222.8, 

1222.5* 

1232,1, 

1231,8* 

1241.3, 

124l.lt 


182.86. 193.21 


1147.5* 

1146. 6 

1157. C* 

1156.2 

1166.3* 

1165.7 

1175.7* 

1175.1 

1185.0* 

U04.5 

1194.3* 

1193*9 

1203.6* 

1203*2 

1212.9* 

1212*5 

1222.2, 

1221.9 

1231.5* 

1231.2 


1240.8, 1240, 6, END 

IFNTHALPY GLYCOL-WATER VS TEMP 
tSAVED TIME STEPS 
ARAD OUTLET TEMPS 

ASECONOARY REGEN COLO INLET TEMPS 

ISECDNOARY BYPASS VALVE POSITION 

SGLYCOL OUTLET TEMPS FROM CONDENSER 

ASFCONOARY REGEN HOT INLET TEMPS 


66. SPACE. 200, END 
67. SPACE. 200, END 
68, SPACE, 200. END 
69, SPACE, 200. END 
70, SPACE, 200, END 
71, SPACE, 200, END 
72. SPACE, 200, END 
73. SPACE. 200. END 
74, SPACE. 200, END 
75, SPACE, 200, END 
76, SPACE, 200, END 
77. SPACE, 200. END 
7ft. SPACE, 200, END 


ASECONDARY REGEN GAINS. SYS 1 
ASECONDARY REGEN GAINS, SYS 2 


ASECONDARV REGEN GAINS, SYS 3 


ASECONDARY REGEN COLD SIDE DELAYS 
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7q.<;PftCe.20C,END 
flC, SPACEt200. EMO 

81,SPACE»20,EMO tSECQNOARV REGEN HOT OUTLET TE«PS 

P2,SPACE»2CtEND 

63t<i*>ACr,20,EN0 

84 t SVSTEH FLOW RATE VERSES RAO PRESS OROP 



6* 

,84. 7. 

10. ,81.5, 

15. ,78. 2, 

20. ,75.0, 

25. ,71.7 

30. t68.4. 

35. 

,65.2, 

40. ,61.9, 

45. ,56. 7, 

50. ,55.4, 

55. ,52.2 

60.r^8.9, 

65* 

,45.6, 

70, ,42.4, 

75. ,39. 1 , 

8C. ,35.9, 

85. , 32.6 


92.,29.3« 95.,26.1tlOO.t22.8»105.,19.6,llO.,16.3,US.a3.0 
120.. 9.a,X2S.f 6.5,130.. 3,3,ENO 
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MAIN PROGRAM LISTING 


130 


110 


Bcn 3EXt'CUTtON 
OOMMON/ALPHA/ WH2n26ni 

CnM«ON/BPAVO/OCONO(31* DMASSm* WTH?30m 

CO»iMnN/CHARLE/J*M,N 

CriHMrtN/DELTA/ IMOEXl 

Cn«M0N/FCH0/ IPQINT 

DIMENSION XdOCOI 

dimension XKU) 

FOtUVALENCE < XK( 1) tX( 1 U 

NOIM=100C 

NTN=0 


SCLDEPUAtKlI iSCAte VISCHSITy 

SPI.irtK2,Al«-t*A%«-t.A6«^ll ASEPARATE DENSITY FROM TEMP 
ARYINV(K2tA6Fll tINVERT DENSITY 

JDIN(K2tA5H,A6 + l,Al*l( $ INVERSE DENSITY VS TEMP 
SPLITIK2, A2*l ,A5*-l, A6FU ^SEPARATE CONDUCTIVITY FROM TEMP 
SPLIT(K2,A3+l,A5*l,ATHI tSEPARATE SPECIFIC HEAT FROM TEMP 
OIVARY(K2t46+l ,A7+1, A7+U ACONDUCTIVITY f SPECIFIC HEAT 
JQIN(K2, A5*l, A7+1, A21 + 1I SCONDUC T I V I T Y/SPEC I F IC HEAT VS TEMP 
SPLI TIX2, A4fl ,A5 + 1, A6+1I ISEPARATF VISCOSITY FROM TEMP 
ARYINVU2.A6+11 S INVERT VISCOSITY 

JQINIK2tA5*l.A6+l*AA+ll SINVERSE VISCOSITY VS TFMP 
SPirTIK9l6»A2A<-l»AA7*UAAS+l) TSEPARATE TEMP FROM ENTHALPY 
J01NIK916,AAfl+l*A47+l,A49+lI TENTHALPY VS TEMP 


REM 

REH DATA ARRAY IS NOW IN PROPER FORM FOR PROGRAM OSAGE 
REM 

REM THE FOLLOWING INITIALISES THE ENTIRE SOLID TEMPERATURE 
REM NETWORK TO COMPENSATE FOR VARYING INITIAL CONOITIONS, THE 


REM TEMPERATURE VALUE IS INPUT AS CONSTANT NUMBER K598e 
STFSQS (K598tK599, TTOl ) 

IPOINT == I 
INDEXl * 0 

J » XKI390) ♦ OeOOOl 
M * XK(399> ♦ OeOOOl 
N = XKUOOl ♦ 0*0001 
00 IOC I = J# Mt N 

WH2026m = (XKU ♦ 325)/KKI! ♦ 3261) - XKI! + 325) 
CONTINUE 


f%eM«A^P«BEGlN INITIALUATION OF SECONDARY LOOP COMPONEMTSA***^ 
REM FIRST INITIALUE CONStANT LOCATIONS 


f 

F 

f 

F 

F 

F 

F 

F 

F 

F 


F 

F 

F 

F 

F 

F 

F 

F 


on no !*JtM,N F 

REM INITIALIZE SEC BP VLV STAT POS TO INITIAL DYNAMIC VALUE 
XKf I^456)=XK(I+374) F 

XKni'365)*XK(H'33l) F 

RTEST=XKU^33ll F 

OlDEGl (RTEST,AlO,TTESt) 

XKIU374)*TTEST F 

XKn^^38l)*U*/XK(U5n F 

REM SET REG HOT INLET ^ CONO GLY OUT 
XKfU39l)*XKn^3e8l F 


REM SET REG COLO INLET = RAO OUTLET 
STFS0S(T4ie»ltK890* 

$TFS0StT636»ltK69l) 

STFS0S(T654tl,K692l 
REM TIME STEP 

STFSQS(.002t<913tA504^l> 

REM NOW SOLVE SYSTEM IN STEADY STATE 
00 30C I=J»H,N F 

lFir*2l I20»l30fl40 F 
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I2r 


no 


no 




147 

150 

155 


160 


165 


170 


200 


03DEG1 (K590,K896.K6f A38*«TEST1 
D!DEGMK390,A3iSTESTI 
OlOEGUK396fA3. TTfiSn 
GO TO 145 

03DEG1 (KS91,K897,Kr ,A38,RTEST I 
D10EGMKS91,43tSTEST) 

[)10EGnK8 97,43,TTeST) 

GO TO 145 

D3DEC1 (K892tK898,K8i A38,aTEST) 
DlDFGl(K892,A3f STESTl 
DlDBGUK8«8,A3tTTEST) 

FLOWC = ( l.-XKU 1-374 n*KK( 1*5 ) 

■JA = 9TFST 

CPC=STEST 

CPHsTTEST 

OC=FLnwC*CPC 

CH=XK( f*5)*CPH 

8ETA»CC/CM 

RNTU=UA/CC 

ETsEXP(-RNTU*(1.-8ETAU 
EFF*t I .-KTl/( 1 ,-BETA*ET) 

P,TFSTsl.-EFF 

STE5T*eFF 

TTE5T*EFF*BETA 

UTESTsl.-TTEST 

XK( I+377)sRTEST*XK( 1+3851 + STEST*XK( I +391 ) 
Xkt t+394|=TTESr*XK( I+385I + UTEST^XM 1+391 | 
IF! XKI 1+3741 .GT.XKI 3811) GO Tn 147 
VTESTsl l,/FLnMC)*67.*XM 361) 

GO TO 150 
VTFST=10C0. 

CONTINUE 

IFtl-2) 155tl60tl65 

STFSQS(RTEST,K9I2,A66+U 

STFS0S(STESTt<9l2,A67+l) 

STFS0SlTTEST,K9l2,A68+n 

STFS0S(UTEST,K912,A69+U 

STFSQS{VTEST,X9l2,A78+l) 

GO TO ITO 

$TFS0S<RTE$T,K412.A70+t) 

STFSQS($TEST,K9l2,A7l+ll 

STFS0StTTFST.K912,A72+l) 

STFSQSIOTEST,K9l2,A73+ll 

STFS0S(VTESTtX912,A79+l) 

GO TO 170 

STFSQSIBTEST,K9 12, 474+1) 
STF$0S(STEST,K912,A75+U 
STFS0S(TTEST,K912,A76+1) 

STF$0SIUTEST,K9 12,477+1) 

STF SOS (VTEST,K9 12,480+1) 

CONTINUE 

RFM INITI4LIZE RENAINOER OF CONSTANTS 
REM CONO GL7 IN * REG COLO OUT 
XKtl+403)*XKtt+377) 

REM INITIALIZE REMAINING ARRAYS 
IFII-2) 200,210,220 

STFSQSIKB90,K914,A51+1) iRAO OUT 

STFSQStK890,K9l2, 454+1) $REC COLO IN 

STFS3SIX879,K912,A57+1 ) *BP VLV FRAC 

STFSQSIK893,K9I5,A60+1 ) tCONO GLV OUT 
$TFS0S(K893,K912,A63+n tRFG HOT IN 

STFSQStK899,K914,ABl+l) iREG HOT OUT 

GO TO 230 


F 


F 


F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 


f 


F 


F 


F 

F 


F 
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2 i: 


GO 


220 


STF 

5QS 

(K891 

.K9U 

t AS2<^ 

11 

STF 

SOS 

(K891 

,K91? 

f A5S+ 

I I 

STFSQS 

1«880 

»K912 

tA50 + 

1 I 

STF 

SOS 

(K89A 

,K9l'> 

t A6l*^ 

l> 

STF 

SOS 

(K894 

fK912 

1 A64^^ 

n 

STF 

SOS 

(K900 

tK9l'f 

tA82 + 

U 

TO ?30 





STF 

SQS 

U892 

»K91^ 

# A53+ 

1 ) 

STF 

SOS 

{K892 

.K912 

t A564^ 

U 

STF 

SOS 

1K881 

,K912 

tA59^ 

n 

STF 

SOS 

(K89S 

.K915 

» AA2 + 

n 

STF 

SOS 

(K895 

.K912 

t A65 + 

n 

STFSOS 

(K901 

.K91<» 

, A83 + 

L I 


230 CONTINUE 

REH END OF I N I T t At I / AT I ON L OOP 
3C0 CONTINUE 

VAR0L2 1 
CNFRWO i 


ANALV2E PRESSURE NETWORK AND OBTAIN FLOW RATES 
PERFORH TRANSIENT ANALYSIS 


END 

BCD ^VARIABLES I 
CnWMON/AtPHA/ HH2026(3I 

CORRON/BRAVO/ClCONOmtOMASS(3ItWTH230I3I 


F 


F 

F 


F 

F 


cqmmon/charle/j»m,n 

COMMON/FOXTRT/ PCKOHm 

DIMENSION SVnLm3ltWH202P(3)fPRES3U3)*SVOL3lO) » 

I WH 2030 m,SV 0 L 21 ( 3 UDSPVnLm.SFCH 2 ( 3 I,SFC 02 ( 3 lt 

? OGENATI 3 »tOH2( 3> t002(3* *0ELECT( 31tPARPH2( 31 tSVOLSaO > , 

3 WOTCP2t3 I ,QRS( 3 U QSTORO 13 J , SVOt 1 2 f 3 ) * QSM( 3> , 

4 QREACT(3)» 0H2O(3U QnUTt3)t QlN(3)t OELTAOI3I* 

5 TSKTNlf 3>tDTEMP(3) »TSKrN2(3It WH2D20(3)t DIFFIBIt 

6 WDTCPU3)v WH2G25I3I 
DIMENSION XKU) 

EQUIVALENCE (XKm,KUM 

REM INTERPOLATE 00 LOOP INDEX ARRAYS 

DlOEGUTfMENt A26,K902) % J 

DIDEGMTTMEN, A27tK903l S M 

DlDEOl (TIMENf A2RfK9041 S N 

REM CHANGE INDICES FROM FLOATING POINT TO INTEGER 

FlX(Kg02»K26> t J 

FIX(K903tK2d> i M 

FIXIK904,K32^ % N 

J»K(26» 

M*K<28I 

N=K(32» 


F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 


f 

F 

F 



DO 


687 I 

- U3 


667 

XK( 

I 

♦ 469) 

» XKI I 

^ 51 


DO 


668 I 

» It 3 


668 

XKI 

1 

+ 5) » 

O.COOl 



on 


669 I 


N 

669 

XKt 

I 

♦ 51 = 

XKtl ♦ 

469 J 


DIOEGI ITIMEN* A4l«K9S4) 
OlDEGUTlMENt A42,K85S) 
01DEGUTIMENtA43tK856> 



DO 

500 I = Jt M 

* N 



DT 

* 10* 




IF 

(1-21 1. ?f 

3 


1 


niOEG2(K336t 

A39* 

K920I 



tflDEr>2(<836t 

A40f 

K932) 


GO 

TO 4 



2 


01DEG2(KS37t 

A39« 

K921) 



010EG2tK83r, 

A40t 

K933I 


GO 

TO 4 




3 OlDEG2tKe38tA39»K922l 

OlOEG2IK83at A40t K934J 


F 

F 

F 

F 

F 

F 

% SYSTEM I 
i SYSTEM 2 
% SYSTEM 3 


F 

F 

F 


F 


F 
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rONTINDF F 

SVOLUt I) = XK( 1*^15) F 

OftftTE = 0.C5 F 

I OLD = I F 

nn 90 jj = It 100 F 

riM202P(n = XK(399)«XKn+349)*DTtMEU F 

MH3n25(n = WH2Q26(t) XK(I + 345 ) *WH202PU ) F 

PCKDHd) = XK( lF325t/( XK( 1+325 ) + HH2025n)) F 

XK( 1+4241 = PCKOHt I) F 

IF (I - 21 21, 22. 23 F 

21 02DEr.2(K929t <847, A44, <9261 

GO TO 24 F 

22 020EG2(K930t <848, A44, <9271 

GO TO 74 F 

23 D2DEG2 (<931,<849,A44,K928l 

24 CONTINUE F 

PRES3K I ) = X<( I+42U F 

SV0L3KI) = (X<( 335t*( X<( I+3421+460. n/(PRES3l( n + l44. » F 

DMASSm = U.O-XKI t + 34511*MH202P( I) F 

WH2030I I I » X<t I+336I+DTIMEU/SVOL3U I 1 F 

WH2020m = WH2030m - DMASSm F 

$VOL2im = X<( I + 336140T INEU/MH2n20( n F 

OSPVQUIl = SVOLlim - SVOL2im F 

IF ( ABSI OSPVOH I 1) - X<( 353)l 100, 100, 35 F 

35 IF (DSPVOLim 45, 100, 50 F 

45 I MULT * I F 

I NEW = I OLD F 

lOLO = -I F 

GO TO 60 F 

50 IMULT * -1 F 

INEM * lOLD F 

into = I F 

60 «(l + 3451 • XKII ♦ 3451 + ORATE*IMULT F 

IF tlNEM - lOLOl 70, 90, 70 F 

70 ORATE = ORATE/2, 0 F 

90 CONTINUE F 

WRITE 16, 951 F 

95 FORMAT ( 79H ITERATIVE 5CHEME DID NOT FIND A SOLtlTlON WITHIN SPECIF 

IFIEO NUMBER Of ITERATIONS! F 

100 CONTINUE F 

WH2026m * WH?"***!!!! 

XKtI+3531 ■ 5.5 ^344 - 1 . ' , ,i 3+ X M , ' ■ ) 

1 +2,654927E-3*Xr.( 1 +349 ) - l.Slj*-. i- ^X<(I+3t ) - 

2 +9, 69469E-2*XK( 1+3421 ^ . 986233fc- i x - ( +349 1 ♦< (1+ 42 i ^ 

3 - 3.573639E-79XKI I + 349I ^<^2+XKU+34? ) - I . 525995c -5 *a<I; +3421 **2 F 

4 -4.8893876-6+XKI I+3491*X<U+3421**2 - L . 292758E- 7«X< (I +342 1 **3 F 

XKII + 3531 * XKII + 3531 - 24.1*1 ,75 - PCKOHdll F 

XKd + 4121 - X<d + 349|*XKd + 3531 F 

SFCH2(I1 * S.292E-5*3t.«XKd+349l F 

SFC02dl * 7,94*SFCH2dl F 

OGENATdl > 51600. *SFCH2d 1 F 

QH2dl » SFCH2I 11*3, 43*1 XKd ♦ 3421 - XKII ♦ 46311 F 

002m = $FC 02 d l*0,2l9*tXKI I + 3421 - XKd ♦ 46611 F 

QELECTdt = 3.41276*XKd+3491*XKd+3531 F 

PARPH2dl » 60. - PRE$3ldl F 

SV0L32(I1 = (XKI3361*IXKd + 3421 + 460. 1 1 /( PARPH2 d 1 *144. 1 F 

WDTCPldl • IXKd+336l/$V0L31dl)*0.445 F 

WDTCP2dl = IXKd+336>/$vaL32d)l*3.48 F 

ORSdl » IHDTCPim+WDTCP2mi*IXKd + 3421-XKIl+3391 1 F 

OSMd) « 3.25*KKd ♦ 3421 - l.6*XKd ♦ 4391 - 882. F 

OSTDRDdl * OGFNATIIl - QH2M1 - Q02II1 - OELFCTdl - QRSdl - F 

1 OSMlIl F 

KKd+3421 = XKd+3421 ♦ tOSTDRO d 1 *OTl MEUl /XKI 358 I F 
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ISC 


160 

US 

168 

170 


10C- 


SCO 


(XK( 336 t*( XK( 1+031 t +660. ) >/( PARPH2( 11*166.) F 

( XK( T+3361/SV0L121 II -XK U +36 9 ) *XK { 663 ) ) *0TIMEU F 

(XX( 665 l*OM4SS( I n/OTIM£l) ♦ F 

XK(6661*(XKU + 369t*XM369) - OM4S S( I I /DT I f*E U 1 F 

s (OM4SSM )*(XK( 1 + 331) - XKtl + 639) )*1.0)/r)Tf "1EU F 

= 3.66*WTH230( I )/nTIMEU + XK( I +336 ) /XK U +6 1 5 ) • 1 . 0 F 

+ DM4SSU l*l.O/DTIMeu F 

= ll60.*D+t4S$n l/DMMEU F 

« n XK(I +336)/SV0L 12(11 )*3.66+( XX(l+336)/SVnill((ll*l.C)/ F 

( (WTH2 30( I )/0TlMElJ)*3.66+(WM2O30( 1 ) /O T I MEU ) * 1 . C ) F 

» -1.2*XK( I + 362 ) + 882. F 

371) * XKU + 362 ) - (XKU + 339) - XK(t + 33 1 ) ) *CONS C 3 F 

> (XK(I ♦ 371 ) - XK(I + 331 ))*C 0 NST 1 + Cr)f 4 ST 2 F 

QREACTm F 

* OCOND(I) + OH 2 nU) ♦ QELECTtn + Os'll I ) ♦ OSrORnU) F 

DELTAO(I) S OOUT(I) - OINMl F 

IF ( 4 BS( 0 ELT 4 QI I ) ) - XKU 69 )) 180 , 180 , 160 F 

F 


svnLi2(i) s 
WTH230m s 
QfiEACni) = 
I 

0H2nn ) 

CONSTl i 
i 

COMST2 ^ 
CDNST3 ’ 

I 

CnNST6 i 
XKU + 

ocoNni 1 1 
oourm = 
ell'll i ) 


180, 180, 160 
IF (4BSrCnNST6 - XK(I + 339)) - XKI650)) 180, 

IF (ABStCGKtSTA - CGNSTS) - ABSIC0NST6 - XK(l ♦ 
DT = -OT/2, 

C0NST5 » XK( I + 339) 

XKU ♦ 339) « XKU ♦ 339) ♦ OT 

GO rn 150 


DM4SSI I ) /DTIMEU 
XK( I+365)*100. 


180, 165 
339) )) 16R,168,170F 

F 

F 

F 

F 

F 

F 

F 

F 


CONTINUE 
XK( 1*650) 

XK( 1+653) 

CONTINUE 

REM INTERPOLATE BYPASS OPTION 

D10EG1(TIMEN,A35,K905) I SYSTEM 

DIOFGKTIHEN, A36,K906) $ SYSTEM 

010E6HTIMEN,A37,K907) » SYSTEM 

PEM CHANGE BYPASS FROM FLOATING POINT TO INTEGER 
FIX(K905,X9) i SYSTEM 

F1X(K906,K10) % SYSTEM 

FIX(K907,K11) $ SYSTEM 

REM INTERPOLATE HEAT FLUX ARRAYS 
RTEST - TIMEM 

IF(TIMEM.GE.1.55) RTEST * TIMEM - 1 . 55*FL0AT ( I F I XI T I MEM/ 1 . 55 ) I 
OllCYL IK600,RTEST,A13,K12) 

OllCYL tX600,RTEST,A16,K13) 

OtlCYL(K600,RTEST,A15,K16) 

OllCYL IK600,RTEST,A16,K IS) 

OllCYL I K600, RTEST, A 17, K 16) 

D11CYLIK600,RTEST,A18,K17) 

011CYL(K600,RTEST,A19,K18) 

OllCYL I K600, RTEST, A20,K 19) 

APPLY INCIDENT HEAT 


ARRAY 


POINT 

1 

2 

3 


PANEL 

PANEL 

PANEL 

PANEL 

PANEL 

PANEL 

panel 

PANEL 


REM 

REM 


ARYHPYU20,A12 + l,Kl2,070t) iPANEL I 
ARYMPY|K20,A12+1,K13,Q725) iPANEL 2 
ARYMPYtX20,Al2+l,Kl6,Q769) APANEL 3 
■ARVMPY(X20,A12+1,K15,Q773) APANEL 6 
ARYMPY(K20,A12+1,K16,0797) IPANEL 5 
ARYMPYtK20,A12+l,Kl7,Q82U IPANEL 6 
ARYMPYtK20,A12+l,K18,QS65) iPANEL 7 
ARYMPVtK20,A12+l,Kl9,a869) IPANEL 8 


REM 

REM 

REM 


BUILO A TEMPERATURE ARRAY TO EVALUATE THERMAL FLUID FLOW COND 

BLDARY(A9+l,T60l,T602,T603,T604,T605,T606,T607,T608,T609 
T6 10, T6U ,1612,1613,1616,1615,1616,7617,1611 
T619,T620,T62l,T622,T623,T626.T625,T626,T627) 
BL0ARYIA9 + 2 8,T628,T629,T6 3i>,T631,T632,T6 33 ,T636,T635,T629 
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d 

d 

d 

d 

d 


d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 


d 
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5 


10 


11 


12 


REM 


REM 

PEm 

REM 

REM 

REM 

REM 

REH 

REM 

REM 

REM 

REM 

REM 


COMPUTE (K-FACTOR )/(2*OENSl rv*t A**2) ) 

DIOIMI (A8, AR+l, Al.KTOl, A23+1) 
niVARY(A23tCS0l,A22*l»G501t 
ADDARY (A23.A23* I, G50l,G901>. 

- ALL CONSTANTS FOR PRESSURE DROP HAVE REEN COMPUTED — SUM 
FOR RESPECTIVE SYSTEMS ANO SOLVE FOR DELR 


IN THE FOLLOWING SYSTEMS COMPRESSIONS A FLAG HAS REFN SET 
TO INDICATE FLOW DIRECTION THROUGH THE BYPASS VALVES — THE 


SYSTEM CONDUCTANCE 
FLAG FOR SYSTEM 

Flag for system 

FLAG FOR SYSTEM 


IS CALCULATED ACCORDINGLY 
1 * KP tO-RAOIATOR t 1-ftYPASS) 
7 » KlOtO-RAOIATOR t 1-8YPASS) 
3 » KU (O-RAOI ATOR , 1-flYPASSt 


COMPRESS PRESSURE NETWORK FOR SYSTEM \ 


SUMAPY(K3*G50l.K582) 

SUMARY(KA,G516tK583» 

IFI KIR) .EQ.O) Kl ll7)aK(116) 
IF(K(9I .EO.DXKI U7t*Clft6) 

AODIK582»K584,G526,K5B5I 


REM 

REM COMPRESS PRESSURE NETWORK FOR SYSTEM 2 
REH 


SUMAR¥IK3,G528»K586I 
SUMARY(KA*G5A3,K387t 
IFtKt IC) .EQ.O) Kt 121 )>Kl 120) 

IFIKIlO.EQ.ltXKI 121)>G( 113) 
A0D(K58G,K58a»G553,K589l 

REM 

REM COMPRESS PRESSURE NETWORK FOR SYSTFM 3 
REM 

SUMARY(K3»G555iK59D) 

SUMARY IK4, 6570» K59 1 ) 

IF(Km).EO*0) K(125)«Ktl24) 

IPIKI 11) .EQ.l )XKt 125)>Gt 140) 
ADD(K590,K592,G580,K593) 
SCALEtK25,KS85*K98StK5B9tK589,K593»K593> 
ARIN0V(3tK6tl.tK821) tINVERT FLOW RATES 
REH SAVE TIME STEP 


REM FIRST TIME STEP» OTIMEU*0.» SO DON'T STORE IT 
IFIDTIMEU.EQ.O. ) GO TO S 
SLPARY(0TIMEU,A50) 

CONTINUE 
bo 700 1>J*M,N 
IFtI-2) 10*11,12 

SLPARY<T6ia,ASl) 

01DEG1IT618, Al.RTESri 
MLTPLV1K821,RTEST,K864,TTEST) 
DELAYl(A50,A5l.TTEST,K890,lTEST) 
SLPARY(iK890,A54) 

ARYST0tK9l2,VTEST,AS7+l) 

GO TO 13 

SLPARYtT636,A52) 

01DEG1IT636,A1,RTEST) 

MLTPLVIKB22,R TEST, K 864, tT£$T) 

DELAY! (A30*A82 ,TTEST,KS91, 1 TEST) 
SLPARYtK891,AS5) 

ARVSTO(K912,VTEST.A59*l) 

GO TO 13 

SLPARV(T684,A53) 

01DE61tT654,Al,RTEST) 

MLTPLY(K923,RTEST,K864,TTESr) 


F 

F 


F 

F 


F 

F 


F 

F 

F 

F 


F 


F 


- 67 - 


SD 70-266 



Space Division 

North American Rockwell 


DELAVl tft50,A53tTTEST.Ke92,!TEST) 

SLPARy A56) 

ARYST0(K912, VTE$T,A59*1) 

13 CONTINUE F 

REN OETERNINE BYPASSED fPACTION 

RTEST-KM I V53U P 

STFST = XK ( 1+365) 

VTEST = XK(I+456) F 

lF(Kn+36B)) 2C,30f3D F 

20 IFtftTEST.LE*5TEST) GO TO 35 F 

25 D1DEG1(RTEST,A10,UTE$T) 

tFtUTEST.LE.VTEST) GO TO 39 F 

K(I+36B)= 1 F 

DIDEGKRTFST.AIO.UTEST) 

VTFSTsUTEST F 

GO TO 39 F 

3C IFIRTEST.GE.STESTl GO TO 25 F 

35 D10EGl(RTEST,All,UTESn 

IFIUTEST.GE. VTEST) GO TO 39 F 

K( I +368)» -I F 

OlOEGltRTEST^AlUUTESTI 

VTEST=l)TEST F 

39 CONTINUE F 

XKU+A56) * VTEST F 

IFM-2) A0rAl,A2 F 

AO SLRARY(VTEST« ASTI 

CNVLTNI A50,A57,0,2,Ka79) 

GO TO A3 F 

A1 SLPARYlVTESTt A5B» 

CNVLTNI A50, A58.0, 2, K880) 

GO TO A3 F 

A2 SLPARYtVTEST»A59l 

CNVLTNI A50,A59,0,2,K88ll 

A3 CONTINUE F 

REN OETERNINE COLO SIDE REGENERATOR FLOW 
FLDRGC * I l.-XKI I+37An*KKt 1+5) F 

REN OFTERNINE REGENERATOR COLO SIDE DELAY 

REN IF BP FRACTION IS LARGE, SET FLORGC TO A SMALL FINITE VALUE 
IFIXKI I+37AI.GT ,XK( 36U) FLORGC^U .-XK I 38 1 > ) AXK I 1 + 5) F 

TTEST-IXKI 1+385) + XK t I +377 ) ) /2 . F 

01D£GHTTEST,A1,HTEST) 

TTEST=XKI305)*U./FLORGC)*RTEST*XK(361) F 

REN SLIDE DELAY INTO ARRAY AND DETERMINE DELAY, COND, OUTLET TO 
REM REGEN HOT INLET 

IFII-2I 70,80,90 F 

7C SLPARYITT6ST,A7fl) 

DtDEGllK893,Al,RTE$T) 

MLTPLYIK82l,RTEST,K867,TTEST) 

TTE$T*TTEST-OTINEU F 

REN DELAYllA50,A60,TTEST,K896,ITEST) DELETED DUE TO INSTABILITY 
CNVLTNI A50,A60, 0,2, K896) 

KKI392I = XKI392) + A60 , F 

XKIA6A) = XKIA60) ♦ (XKI3921 - XK 1 460 ) ) *< . 53 ♦ .C02*XK(6) ♦ .018* F 
1 XKIA621I F 

DH = IXKIA60) - XKIA6AI )*(.152 * ,000022 5*1 XK I A6A) ♦ XK(A60))) F 

I *XKtA62)/XK(6) F 

X * 116AA0.5 + 2A01.36A*OH + XKt 392 )♦( 682 ,4676 + XKI392H F 

TOW = SQRTtX) - 341.2338 F 

XKIA67I * XKtA61) ♦ .75*(TQH - XK(A6l)) F 

OH * tXKIA61) - XKlA67l)*3,AT*XKIA63)/XKt6» F 

X = 116440,5 + 240l,364*0H ♦ TOW*t 682 ,4676 + TOW) F 

XKI392) * SQRTIX) - 341.2338 - 460. F 

XKI464) = XKI464) - 460. F 

XKI467) = XKI467) - 460. F 
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SLPARY(K8<)6t A63) 

r,n TO 100 F 

80 SLPARV(TTEST,A79) 

niOEGUK894tAl,RTESn 

MLTPLY(K822,RTEST,KR67,TTESTI 

TTEST = TTEST-OTtMELI F 

REM OFLAYUA50,A61,TTEST,K897,ITEST) OECETFO DUE TO IMSTARILITY 
CNYLTNt A'SOt A6lt 0, 2»K897I 

XK(393t = XK(393I ♦ ^60. F 

XKtA65) = XK(460) ♦ tXK(393) - XK K60 ) ) ♦ ( . SB f .002*XK(7) ♦ .OIR* F 
1 XX(462n F 

OH = (XK(A60) - XK(465M*(.152 <■ . 000022 5* ( XK ( A65 ) + XKtAftOlU F 
1 *XKU62)/XK(7) F 

X > 116440.5 * 2401.364*OH + XK ( 393 )*( 682 .4676 + XK(393)» F 

TOM = SQRT(X) - 341,2338 F 

XK(468J * XK(461I + .75*(TOM - XK(46in F 

OH > (XK(46U - XK(468) l*3.47*XK(463) /XK( T» F 

X = 116440.5 + 2401.364*DH ♦ TOM*( 682.4676 * TOM) F 

XK(393( s SQRT(X) - 341,2338 - 460. F 

XK(465) = XM4651 - 460. F 

XKt468t = XK(468) - 460. F 

SLPARY(K897, A641 

SO TO 100 F 

90 SLPARY(TTEST,A80) 

D1DFGI(KB95, Al.RTESTt 
*3LrPLVCK8 23,RTEST,K867,TTEST) 

TTEST=TTEST-0T1HEU ' F 

REM OEtAYl(A5O,A62,TTEST,K890,ITEST) OELETEO OUF TO INSTABILITY 
CNVLTNI A50, A6 2.0,2, K898) 

XKt394) a XKI394) ♦ 460. F 

XKI4661 * XK(460( * (XKI394) - XK ( 460 ) ) * ( . 53 ♦ ,002*XK(8) + .018* F 
1 XKI462M F 

OH a IXK14601 - XK(466n*I.lS2 * ,0000225*1 XK 1 466 ) + XK(460m F 
I *XK(462)/XK(8I F 

X = 116440.5 * 2401.364*OH * XKI 394)*( 682.4676 + XK t 394 )) F 

TOM a SQRT(X) - 341.2338 F 

XK(469» a XK(46U ♦ .75*(T0M - XK(46l)) F 

OH a (XK(461t - XK(469)1*3.47*XK(463)/XK(81 F 

X a 116440.5 ♦ 240l.364*OH ♦ T0M*I682.4676 * TOM) F 

XKI394) a SORT(X) - 341.2338 - 460. F 

XKI466) a XKI466) - 460. F 

XKI469) a XKI469) - 460. F 

SLPARY(K898,A65) 

100 CONTINUE F 

REM COMPUTE REGEN HOT SIDE DELAY 

TTEST*(XKU*394» ♦ XK(I*39U)/2. F 

blbEGlITTEST,Al,RTEST) 

XK(33)=XK(385)*XK( I*3l6)*RTEST*XKt 364) F 

PE)^ FIND TERMINAL ENTHALPIES AND UA OF EACH REGENERATOR 
IF(I-2| 110,120,130 F 

110 030EG1(K890,K896,K6,A38,RTE$T) SUA 

_ DlDEGl (K896, A24,STESn tHMl 

b"i0EGliK899,A24,TTE5T) $HH2 

010EG1(K890,A24,UTEST) »HC I 

010eGltX«a2,A24,VTEST) tHC2 

GO TO 140 F 

120 03DEGltK891,Ka97,K7,A38,RTESn iUA 

01DEGUK897,A24,STESTI tHHl 

6ibEGitK9bO,A24,TTESTl *HH2 

01DEG1(K891,A24,UTEST) tHCl 

010EG1IK883,A24,VTEST» 1HC2 

GO TO 140 F 


- 69 - 


SD 70-266 


f 



Space Division 

North Arnerican Rockwell 


13 ? 


UO 


150 


160 


17C 


030rr.l (K892*K898,K8tA3a,RTESTI tUA 

niDFGl (K89B, A24, STEST) $HH 1 

010FGl(K90l,A24tTTFST) tHH2 

DlOEGl(K992,A24»UTESn »HCl 

0l0EGl(Kftft4,A24,VTESTl iHr,2 

CONTINUE 


ftFM COMPUTE TEMPERATURE AND ENTHALPY DROPS ON EACH SIDE 
UAsPTEST 
HHI=STEST 
HH2=TTFST 
HCt=UTFST 
HC?=VTFST 
DHH=AnS»HHl-HH2l 
nHC=ARS(HCl-HC2l 
OTH = ABS(XK( I+OPl l-XKI 1+394) 1 
DTC = ABS(XK( 1 + 385) -XKU +377) ) 

CC=FLORGC«( l./OTC)*OHC 
CH=XKt I+5)*tl./OTH)»OHH 
BETAsCC/CH 
RNTU=UA/CC 

FT=exP( -RNTU*( l.-BFTA) ) 

FFF = ( l.-FT)/( t.-9ETA*ET) 

REM COMPUTE REGEN GAINS BASED ON REGEN INLET CONDITIONS 


RTEST=U-EFF 

STEST=EFF 

TTE$T=EFF*BETA 


UTEST=1.-TTEST 

»EM SLIDE GAIN CONSTANTS INTO ARRAYS AND 
PEM DETERMINE TRANSIENT OUTLET TEMPS 
IFU-2) I50tl60,170 

SLPARY(RTEST,A66) SSYSTEM I 
SLPARY(STEST,A67) 

SLPARY(TTEST,A6R) 

SLPARY(UTEST,A69) 


DELAY2(A50,A54,A78» 
CNVLTE(A50,A54*A66, 
CNVLTEIA50*A63.A67, 
0ELAYl(A50»A63,K.33t 
CNVLTEU50.A54tA68, 
CNVLTE(A50,A63t A69, 
VTFSTsTTEST+UTEST 

SLPAPYIVTEST, A81) 
GO TO IRO 


VTESTtJTESTl 


JTEST, l.RTEST) 

SOUTPUT 

Til 

C,l,STESTI 

JOUTPUT 

T12 

VTEST, JTEST) 



0,l,TTEST) 

SOUTPUT 

T21 

JTEST, 1, UTEST) 

iOUTPUT 

T22 


SLPARYIRTEST.A70) 
$LPARY(STEST,A7l) 
SLPARYITTEST.A72) 
SLPARY(UTFST,A73I 
DELAY2I A50* A55*A79tVTE 
CNVLTEU50,A55.A7O,JTE 
CNVLTEIA50.A64t A71,0. 1 
DELAYI(A50»A64,K33,VTE 
CNVLTEI A50. A55,A72tCt 1 
CNVLTF(A50»A64,A73,JTE 
VTESTsTTEST + UTEST 

SLPARYtVTEST,A82) 


STtJTEST) 


5T,l,RTEST) 

iOUTPUT 

TU 

,STEST) 

iOUTPUT 

T12 

ST, JTEST) 



,TTEST) 

iOUTPUT 

T2l 

ST, ItUTEST) 

♦OUTPUT 

T22 


GO TO 180 


SLPARY(RTEST,A74) SSYSTEM 3 
SLPARYISTESTt A75) 

SLPARY(TTFST,A76) 

SLPARYIUTEST, ATT) 

DEL AY2( A50»A56.A80.VTESTt JTEST) 

CNVLTEIA50,A56,A74t JTEST.URTEST) AOUTPUT Til 

CNVLTEI A50,A65» A7 5,0,1, STEST I iOUTPUT T12 

DELAYll A5D,A65,K33,VTEST, JTESn 


F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 


F 


F 

F 


F 

F 
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CNVLTHt A50, A5A, 176,C, It TTfSrt SflLITPUT T21 

CNVLTEt A50,A65, 477 tJ TEST, l.UTEST ) tOUTPiJT T22 
VTESTsTTEST + UTEST F 

SLPAPVtVTEST, A83I 

REM COLD ♦ HOT TRANSIENT OUTLET TEMPS 
180 XK( H-377) = RTFST + STFST F 

XKn»-394) = TTEST + UTEST F 

TTEST = (XK(I + 39U * XK(f ♦ 3<?4)l/2. F 

niDEGl(TTEST,A3,UTESTt $ SEC REG HOT CP 

TTEST = tXKIT + 377) + XKU ♦ 385H/2. F 

D10EGHTTEST,A3,VTESTt $ SEC REG COLO CP 

XXH + 377) = XKU + 385) * UXKU * 391) - XKI I + 394 ) ) *UTE ST ) / F 

I ttl. - XKU * 374)t*VTEST) F 

IF (XKU + 377».GT.XK( [ + 391)) XKU ♦ 377) = XKU ♦ 39U F 

XKU + 377) = (XK(I ♦ 377) + XKU + 322))/2. F 

XK( I + 322) = XKU f 377 ) F 

RFM COMPUTE DELAY TO RAD INLET AND SET RAO INLET TEMP 
lFU-2) 190,200,210 F 

190 DlDEGUKe99, A URTEST) 

MLTPLYtK82I ,RTEST,K969, TTEST) 

DFLAYl ( A50,A8l,TTEST,T601 , 1 TEST) 

GO TO 220 F 

200 010EG1(K900, A URTEST) 

MLTPLY(K822,RTFST,K869, TTEST) 

OELAVU A50,A82,TT6ST,T619, ITEST) 

GO TO 220 F 

210 OIDEGI(K901,A1,RTEST) 

HLTPtV(K823,RTEST,K869,TTEST) 

0ELAVUA50,A83,TTEST,T837,ITEST) 

220 CONTINUE F 

REH COMPUTE CONDENSER INLET TEMP 

IF( KK< H-374) ,GT,XK( 381) ) GO TO 290 F 

RTEST3XKU4385) F 

STEST=XKU4377) F 


OIDEGURTEST, A24,TTEST) SENTHALPY BYPASS 

OlOEGl (STEST,A24, UTEST) »REGEN ENTHALPY 

260 VTEST=TTFST*XKU4374) + UTESTAU.-XKU + 374) ) F 

0l0EGl(VTEST,A49,RrEST) «TEMP OUT OF VALVE 


XKI r+403>*RTESr F 

GO TO 300 F 

290 XKU+403) = XKU«-385I F 

300 CONTINUE F 

IF UNOEXtl 450, 450, 301 F 

301 IN0EX2 i F 

0TFMP«XK(447) F 

ISGN=.Kt448) F 

XKI I 4 3)65) 3 XKU * 331 ) F 

310 T20LD»XKt I+331 ) F 

T4nL0>XKU43S8) F 

iFU-2j 320,330,340 F 

320 01DEGl(K908,A24,RTEST) 

01DEG1(K893,A24,STEST) 

GO TO "350 F 

330 DlDEGltK909,A24,RTEST) 

DlDEGltJt894,A24,$TEST) 

" GO to “350 ■ ■ ^ p 

340 01DEGl(K910,A24,RTEST) 

OlDE6lIK895,A24,STEST) 

350 CONTINUE F 

IF 1INDEX2 - 100) 354, 354, 445 F 

354 OGLY = XKU + 51*(STEST - RTEST) F 

bfGAS-'XKI 1*371 )-XKI 1*331) F 


351 OGAS » DTGAS*(3,44*MTH230U)/DTIMFU * XK U * 336)/XKU ♦ 415)*1.0)F 
1 *(0MASSU )*ll60,*DMASSt I )*0TGAS* I .0 ) /OT I MEU F 
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XX(T ♦ 331) = XK(I ¥ 388) - XK U t 3 1 8 ) * ( ?2« T')-i9 + F 

I C.687969198*(JGAS ♦ 9 . 0 1 45E- 7’OQC AS **2 I F 

IF tARS(QGLY-OGAS) .LT.XK( A07) I GH TO A50 F 

ISGNSV'ISGN F 

IFOGLV.LT.OGAS) ISGNs + l F 

IF(QGLY.GT.OGAS) ISGN*-1 F 

TINCR"OreMP*ISGN F 

IFI I SGN.eO. ISGNSV) GO TH A40 F 

XKU+33l)=T20LG F 

XK ( I+38S)=TAOLG F 

OTEMP'OTEHP/2,0 F 

INPexZ s INOFX2 ♦ 1 F 

on TG 31G F 

4AC, XKd ♦■383)sXKn 4.383 )+TINCR F 

INDEX2 « INPEX2 ♦ 1 F 

GO TO 310 F 

AAS WRITE <8, 848) F 

388 FOP*XAT UHOt 7AHCONDENSE8 HFAT RALANCE NOT COHPlFTEO WITHIN SPECIFF 
IIED NUMBER OF ITERATIONS! F 

850 OCONOm = OGAS F 

IFII-2) 851*852*853 F 

851 SLPARY(Kfl93,a80) 

GO TO 7G0 F 

852 SLPARY(K898,A81) 

GO TO 700 F 

853 SLPARY(K895*A82) 
lif END OF L 

700 CONTINUE F 

INOEXl si F 

R£M*****A*A8P**8«88*«P8P ****** 8 88** 8*88 8*8 *************** *8 ****** 

FNO 

flCO 3nuT»UT CALLS 

COMMON/RRAVO/OCONOI 3) tOMASSI 3) ,WTH230( 3) F 

COMMON/CHARLF/ J* M, N F 

CDMMON/ECHO/ IPOINT F 

COMMON/FnXTRT/ PCKOHO) F 

DIMENSION ATEMP(50> jl F 

DIMENSION XXm F 

DIMENSION CRTARl(?C0). CRTAR2I200)* CRT AR 3 ( 200) «. CR TAR8 1 200 ) * F 

I CRTAR5I72I* CHTAR6I72), CRTAR7I63) F 

DIMENSION DATEI2) F 

OATA ITlT/0/ F 

DATA CRTAR5/ 'MISS', MON • * ' T I M£ • , ' - •*• MIN'* 13*' '/ F 

DATA CRTAR5/ 'SYST'.'EM T • * ' EMPE ' * ' R A TU* * ' RE S ',' . '*'OEG '* F 
I 'F '♦ 10*' •/ F 

EQUIVALENCE IXKUMKim F 

CALL CAMRAVI9) F 

TIME » 0.0 F 

DO 100 N1 * 1, 63 F 

CRTAR7INI) * TIME F 

TIME • TIME 8 3.0 F 

lOD CONTINUE F 


REM SOLVE FOP RADIATOR PRESSURE DROP 

MLTPLYIK5.K6,K5R5,K8S2) tRADlATOft PRESSURE DROP FOR SYSTEM I 

MLTPLY(K7,k 7*<589,K6»3) tRADIATOR PRESSURE DROP FOR SYSTEM 2 

MLTPLY(XR,K3.K593,K688I SRADIATOR PRESSURE DROP FDR SYSTEM 3 

01DEG1U832,AS8*K6> % FLOW RATE • SYSTEM 1 

I)10Er.llX683«AR8,K7l $ FLOW RATE - SYSTEM 2 

01DEG1(K668,A38,K8) t FLOW RATE > SYSTEM 3 

REM FLUID HEAT LOSS FOR SYSTEM I 
OMTRI 1 16, A9*l ,C80ltA25*l) 

RFM FLUID HEAT LOSS FOR SYSTEM 2 
QMTRI I 16,A9*19*G819,A25*17I 
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REW FLU!0 HEAT LOSS FOR STfiTEM 3 
Q^TR[ U6,49 + 37,G637,A2^>+33) 
[F(Kf9),FQ.n Gn TO 1 

SUGARY U 6 fA 25 *^ltK 80 U 
GO TO 2 

1 CONTINUE 

SUMARY no f 425^^1 *Kd 0 n 

2 Tf(Kno),Eg,u GH rn 3 

SimARV ( 16 1 425^^1 7tKG02) 

GO TO 4 

3 CONTINUE 

SUMARYUC t A 25 <^l 7 tK 0 O 21 

4 IF(KflU*EQ,n GO TO 5 


F 

F 

F 

F 

F 

F 

F 


SUMARYtUt A 25 + 33 ,K 803 I 
GO TO 6 
CONTINUE 

SUMARYI ID » 425^^33 fKROB) 


onWER 

OCONO 


AMPS 

L8/HR 


CURRENT 

rf*C0N0 


VOf-TAGE 

WS-RATE, 


TSt 

3 )(t 


TSE 


VOLTS 

PERCENT# 


DEGF 


DEGF 


OEGF. 


OEGF 

/) 


6 CONTINUE 
CALL DATOUT 
WR!TEf6#m 

U F0RMAT<5)(#99HSV$TFM 

1 TCIP TCEP 

2 12HT02 TH2I 

WR!TEI6f 12) 

12 F0RMAT(l4Xt93HWATTS 

1 DEGF 0TU/HR L8/HR PERCENT# 3X# 12HOEGF 

on 13 i=j#M#N 

13 MR I TE 1 6# 14) I#XK(n4l2»#XKtn^349)#KKU + 353l#)Oai+339),XK(U342I# 

1 XKm 37 l),XK(Iv 33 l)#OCnNDm#XKm 450 )#XKtn 453 )# 

2 XKU ^ 4631# XKir ♦ 466), PCKOHI H 

14 F0RMATI6X«I2#4X#2F8.2, 1X#5F 2 #F 10. 2 , F 1 0 • 3# 2F9. 2 #F8 . 2 #F9. 3 ) 

WRITE (6# 15) 

15 FORMAT I//) 

WRlTEt6;i6) 

16 F0RMATt5X#ll4HSYSTEM 
ID TSRCr T$RCE 

WRTTE(6#m 

17 FORMAT (16X,95HLB/HR BTU/HR OEGF 

1 OEGF DEGF OEGF DEGF 

00 18 iO#M#N 

18 WRITE (6# 191 I#XKU^5)# XKU+296),TI I+243) #T( l>240) #XK( H'212) , 

1 XK( I*3 85), XK( U377>#XK( I +39 U # XK 1 1 ^394 1 #XK(I+403)# 

2 XK(t + 385),KK(I+374» 

19 F0RHAT(6K#r2#5XtF8.2#lX#2F9.2#Fl0*2#F6*3,6F9.2#F8.4) 


FLOW RATE 9RA0 
TSRHI TSRHE 


TRAOfN 
TCIS 

DEGF 
OEGF# /) 


TRADDUT 


F 
F 

F 
F 
F 
F 
F 
F 
F 
F 
F 
F 
F 
F 
F 
F 
F 
F 
F 

DPRAF 


TCES 


PSI A 


8PFS) 


OEGF 


F 

F 

F 

F 

F 

F 

F 

F 

F 



WRITE ( 6 . 20 ) 





F 

20 

FORMAT t //Tit, izOHPANEL 

1 

PANEL 2 

PANEL 

3 

PAF 


INEL A PANEL 5 


PANEL 6 

PANEL 7 


PANEF 


2 L 9 . /» 5 X» 123 MIN 

OUT 

IN OUT 

IN 

OUT 

TF 


3 N nut IN 

OUT 

IN OUT 

(N 

our 

IN F 


A OUT , / » 





F 


DO 21 n « I 93 t 208 





F 


ji ■* 2e<) - li 





F 

21 

ATEMPUl) « Till) 





F 


on 22 12 • 209 , 22 A 





F 


J 2 * 2 A 1 - 12 





F 

22 

ATEMPU 2 ) *> Tf 12 ) 





F 


DO 23 13 * 225 , 2 A 0 





F 


' J 3 i " ‘273 13 





F 

23 

ATEMPIJ 3 ) » TtI 3 ) 





F 


WRITE ( 6 # 25 ) (ATEMP(KT)# XT 

- 1 , 16 ) 



F 


25 


WRITE (6# 
WRITE (6# 
FORMAT (IX# 
NP * 125 


251 lATEMPtKT)# XT - 17# 32> 
25) UTEMP(KT)# XT » 33# 48) 
16F8.1) 
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CRTARl f IPOrMT) = XK(3A3) f- 

CRTARZUPOTNT) = XKn32) F 

CRTARSnPOTNT) = T(24A» F 

CRTAPAUPniNT) = T(2A-lt F 

IF I IPOINT.eO.NP) 00 TO 30 F 

IPOINT = IPnii^T + I F 

GO TO 50 F 

30 CONTINUF F 

IF utit.eq.o) call CDATEVIDATEI F 

IF ( IT[T.GT*0) go to 5D F 

CALL GRF3TV(-1 »3etCRTAR5,CRTAR6,-ft3,CRTAR TtCRTARU 1 ) »2C0, iC, , F 

1 500. ,C. I F 

CALL RITF?VI331,lOlO,l023,ROt2.2C,lt*SYSTE'1 I - 30 AMPS'tl^LY) F 

call PRIMTV(S.DATEt9^3,P95» F 

CALL G9F3TVI 0 » 63 , CRT AR5 1 CR T AR6, -63 ,CRT AR 7 »CR T AR 2( U , 200 . , C . * F 

1 500. fO.) F 

CALL GRF3TVI 0 , 55 , CRT AR5 , CR T AR6, -63 t CRT AR 7.CR T AR 3 ( 1) » 2C0. . 0 . . F 

1 500. ,0.1 F 

call GRF3TV( 0,4A,CRTAR5,CRTAR6,-63,CRTAR7,CRTARA( U ,200. ,C. , F 

1 500. ,0.1 F 

CALL GRF3TV(-l. 3ft,CRTAR5,CRTAR6,-63,CRTAR7,CPTARl(63J ,200. .0, , F 

1 50C.,0.» F 

CALL RITE2VI 331, 10 10, 1023,90,2,20, 1, * SYSTEM 1 - 25 AMR<;<,IrlYI F 

CALL PRTNTVIS, DATE, 963,995) F 

CALL GRF3TVI 0 , 63 ,CRTAR5 , CR TAR6 ,-63 ,CRT AR 7,CR T AR2 ( 6 3) ,200.,0. , F 

I 500. ,0.1 F 

CALL GRF3TV( 0 , 55, CR T AR5 , CR TAR6 , -63 , CRT AR 7 ,CRTAR3 ( 6 3 ) ,200. ,0. , F 

1 500. ,0.) F 

CALL GRF3TVI 0 , 46, CRT AR5 , CRT AR6 , -63 , CRT AR 7 ,CR TARA ( 63 ) ,200 , ,0 . , F 

1 500. ,0.) F 

ITIT » 2 F 

50 CONTINUE F 

END 
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USER 


0001 

0002 

0003 

000^ 


0005 

OOOfc 

0007 

000.9 

0009 

0010 
00 n 
0012 
0013 
OOU 

0015 

0016 

0017 

0018 
OOlt* 
0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 
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SUBROUTINE LISTINGS 


SUBSmjT I-MF CNVLTNfDT , VI ,nST,NFCN,V2 1 

0GU3LF PPFC ISrON OTOP , OTNS . VI J , V IJ PI , V4 V , TMT , T^' rn . V2UP t x $F F f-: , 
1 SCaLF,RrSfl,PMXSB 

OIMENS ION on 1 1 , Vl( U 
E0UIV4LENCE{0,N) 

C 

C VI MUST RE DIMFNSrnNEn NOT LONBEK THAN 

C (DIMENSION OF DTm 

C 

0 = DTU I 
IDT = N 
0 = VI ( 1 ) 

IVl = N 

NT = lOT+l-lVl 
fF(NT) 169, 5t5 
5 IF(NST) 15,15,25 
15 NST = IDT+1 
25 V2DP = 0.0000 
TMT = 0,0000 

MX$S = 10 ‘ 

RMXS& = MX$9 
NSTMl := NST-l 
DO ^0 r=^l,N$TMl 
J = NST-H^l 
□TOP = DT( J) 

JS = J-NT 
VIJ == VUJS> 

VlJPl = Vl( JS + U 
OTNS = OTDP/RMXS^ 

TMTO = TMT 
bb" 70 I S0-l*MXS3 
IF( ISB-n 30f 30f AO 
30 TMT s TMT4^0.5000^0TMS 
GO TO A5 

AO TMT = TMH-OTNS 
A^ RIS9 - ISB 

SCALE • ftlSe/PMXSB 
VAV * VIJPI<^(V1JP1-V1J1*5CALE 
70 V2DP » V2DP^XSFEP(NFCNtTHT)»VAV»DTNS 
J[HT * TMTO^DTOP 
ab CONtlNUE ' ^ ^ 

V2 = V2DP 
RETURN 

\69 WRrTEf6,170) 

170 F0RMAT(A3H error IN DIMENSIONING OF ARRAYS FOR CNVLTN) 

RETURN 
■ END 
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OOQl 

0002 

0003 

ooc^ 


SUBROUTINE CNVL T E ti)T , VI 1 N , NST »NFCN t V 2 ) 

[)OURtE PP EC I SI ON OTDPtDTNStVlJtVlJOl , Vi V * r>J , GJ P 1 tCAV , T VT , T 
1 V2nPtXSFEP fSCiLFf P 1 SB» 

DIMENSION Dim ,V1 m tOAlNt n 
FOUIVALFNCE(DiN> 

C 

C VI AND GAIN MUST BE OIMENSIGNEn EQUAL LENGTH AND MAY 

C NOT 9E longer THAN (DIMENSION OF DTJ+l 

C 


0005 


D = 0T{ 1 } 

0006 


lOT = N 

0007 


D = Vim 

0008 


IVI ^ N 

0009 


D GAIN( U 

0010 


IG = N 

0011 


NT = IDT^l-IVl 

0012 


iFtNT) 169, 5t5 

0013 

5 

NT - IDT<^1-IG 

0014 


IF(NT) I69,10tl0 

0015 

10 

IFHC-tVl) 169,1=;, 169 

0016 

15 

IF(N$T) 20,20,25 

0017 

20 

NST = TDT+1 

0018 

25 

V2DP = 0*0000 


OOIQ TMT = 0.0000 


0020 

0021 

0022 

0023 

002A 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 


»^XSB = 10 

Q ««**««***««**«*«**«** **i»««**<c *««*««*«*«<[«««« *1(1 

RMXSB = MXSB 
NSTMl = NST-X 
nn 80 1=1, NSTMl 
j = NST-I+1 
OTOP = DT(J) 

JS = J-NT 
VIJ = Vl(JS) 

VUPl = VUJS+1) 

6J = GAINUS) 

GJPl = GAIN(JS+n 
DTNS = OTOP/RMXSfl 
TMTO * TMT 
00 70 ISB=l,MXSB 
IF(lSB-lt 30,30,40 
30 TMT = TMT+0.5000*0TNS 
GO TO 45 

40 TMT = TMT+OTNS 
45 RISB = rSB 

SCALF = RISfl/PMXSB 
VAV = VlJPl4(VlJPl-VlJ)*SCALE 
GAV = GJP1+(GJP1-GJ»*SCALE 
70 V2DP = V20P4XSFER(NFCN,TMTI*VAV*GAV*0TNS 
TMT = TMTO+DTDP 
80 CONTINUE 
V2 = V20P 
RETURN 

169 MRITF(6,170) 

170 F0RMAT(43H ERROR IN DIMENSIONING OF ARRAYS FOR CNVITE) 
RETURN 

END 
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0001 

0002 

0003 

000 <* 

0005 

0006 

0007 

0008 
0000 
0010 

0011 

0012 
0013 
OOU 

0015 

0016 
0017 
0013 

0019 

0020 


0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 
0001 
0002 
0003 


0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 
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0016 
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0020 
0021 
0022 
0023 
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SUBROUTINE DELAYKT4R, ARR, TAJ, VAL * INOX) 
niMENSION TARin, ARR(l) 

EQUiVALENCetOtN) 

n ' TARm 

ITAR = N + 1 

0 = ARP ( 1) 

lARR = N + 1 

JTAR = ITAP 

JARR = lARR 

OT = TAU - TARIJTAR) 

10 JTAR = JTAR - I 

IF (JTAR - 2) 100, 20, 20 

20 JARR = JARR - I 

IF (JARR - 2) 200, 30, 30 

30 IF’ (OT) 6'oi 60, 40 
40 OT = DT - TAR(JTAR) 

GO TO 10 
50 INDX = JARR 

VAL = ARR(JARR) 

_R^TURN 

C" 

C ERROR MESSAGES 

C 

100 MRITE (6, 1011 

101 FORMAT (50H OELAYl ERROR--TAU EXCEEDS SUM OF SAVED T[ME STEPS) 
INDX = 2 

VAi. = ARR(2)" 

RETURN 

200 WRITE (6, 201) 

201 FORMAT (28H DELAYl ERROR — ARR TOO SMALL) 

INDX = 2 

_VAL * ARR12) 

WRITE (6," 250) OTIMEU 
250 FORMAT I 12HTHE TIME IS , F10.5) 

RETURN 

END 

SUBROUTINE 0ELAY2(TAR, ARR, TAU, VAL, INDX) 

J^LMENSJON T«IJ^)£^ ARRtll, TAU(l) 

F'QUi valence (D,N) 

C 

C TAR MUST BE DIMENSIONED ONE LESS THAN ARR ANO TAU 

C 

D » TARI 1) 

__ _ITAjR_= N_ 

D"a ARRI 1 ) 

I ARR = N 
0 « TAUI 1) 

I TAU = N 
SUMDT = 0.0 
TAUITG * 0,0 

j' = ITAR + 1 - lARR " ' 

IF(J ) 69,10,69 

10 IFIIARR - ITAOI 69,20,69 
20 bn 5<j I = I, ITAR 
M = ITAR +2-1 
SUMDT = SUMDT + TARIM) 

TAUAV’l = (TAUIM) ♦“TAUIM + l))/2,0 
TAUITG = TAUITG ♦ TAUAVI *TAR ( Ml 
TAUAVC = ( 1,/SUM0T)*TAUITG 
IF I SUMDT - TAUAVC) 50, 60, 60 
50 CONTINUE 

WRITE (6, 51) 
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C0?4 

51 

FORMAT U3H MEAN TAU 

OUT OF RANGE 

HF 

TAR SUM IN 

Of;L AYR ) 

002*^ 

60 

VAl, = JRR(M) 





0026 


INDX = M 





0027 


RETURN 





002S 

69 

WRTTF (6, 701 





0020 

TO 

FORMAT (A3H ERROR IN 

DIMENSIONING 

OF 

ARRAYS FOR 

DELAY?) 

0020 


RETURN 





0031 


END 
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0004 
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0009 
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0001 

0002 

66o3 

0004 
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0000 

0009 
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SUBROUTINE SL *>ARY ( AR YNt ARY 1 

OIWENSTON ARYfl ) 

EOUtVALENCEfOfN} 

O^ARYU) 

IC^N 

00 16 1=2, !C 
10 ARY( [loARYt [4U 

ARY(TC*l >=ARYN 

RETURN 

END 


SUBROUTINE ARYSTOf N , X , A ) 
OlMENStON Aril 
X = A(N) 

RETURN 

END 


nouBLE PRECISION FUNCTION XSFER I NFCN * T I 
DOUBLE PRECISION T 
x'SFER '= T*Tf>do 

C NFCN=1 IS SECONDARY REGENERATOR FUNCTION 

IF INFCN.NE* 1) GO TO 16 
IF(T*GT, 0*52770001 RETURN 

XSFER = I1.000/0*52T70-02>*OEXPI-T/0*52770-02J 

"return 

10 CONTINUE 

C 4r***!0r*tt**^%«4i«««V*««««****«**'««***«*4i««4«*4t*#«4r««!tt*«*«4r««*^«****^* 

C NFCN*2 IS SECONDARY 6 YRASS VALVE FUNCTION 

I F ( ilFCN . N^ 2 1 GO T O 20 
rFU*GT7C,2'e33DOO)" RETURN 

XSFER * a*000/0*29330-02M‘DEXP<-T/0.2B33D-02i 
RETURN 

20 CONTINUE 
RETURN 

END ■ ■ ■ ■ ' ■ 


SUBROUTINE DATOUT 

CQMNON/JF I X CON/ N 

OTMENSinN’NIll 

IF fNf2B>*GE.5B*OR*N<29l*EO*Ot CALL HEAONG 
NI29I = NI2S> +4 

HRifE 16, 1061 N(l>, HiZU N(35i, NIITI, N(36l, NI3T1, 

1 N(30) 

100 FORMAT f/lX t 10 K» * * 4H TTHE, IPE13*5, IX, 6H0TIMEU, 

1"TpF13.5, ‘ix, 7HCSGNINC, 15, IHI, IX I RE 5 , IX , fHOTHPCC (' t 15^ 

2 IH), IX, IPE13*5, IX, 7HARLXCCI, 15, 1K» , IX, 1RE13*5 //) 

RETURN 

END "■ 


SUBROUTINE HEAONG 

COWH ON/F I XCON/ N „ _ _ 

COMHON/TITLE/ H ' ' 

OlMENStON Nil), K(20I 

N(_28l = 15 

Nt291 = NI291 ♦ I 
WRITE (6, 1001 N(29I, H 

100 FORMAT tlHl, 116X, 4HPAGE, 2X, 17 / 4X, 121NSYSTEMS IMPROVED NUHBR 
IICAL DITTFRENCING ANALYZER » STROT" NORTH’^AHERTC^N ROCKW 

2LL CORPORATION - SPACE OtVtSlON, //, 5X, 20A6//I 
RETURN 

■ END 
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0001 



SUeftOllTINE r.RF3TVfLt I ftCD X , BCOV , NP , X t Y t KMA X t XM I N t YMAX tVW I N 1 

GPFOOOlO 


C 



GRF00020 


c 


GRAF3V' PLOTS 3 GRIDS PEP PAGE W THE SC 4020 

GPF0003D 


c 



GRF00040 


c 


444 NOTATION 44« 

GRF00050 


c 



GRF00060 


‘ c 


L NO* OF GRAPHS PER FRAM EU i 2 1 0R3 » * 

GRFOOQTO 


c 


IF L IS NEGATIVE - FRAME WILL BE ADVANCED AND NP POINTS 

GRF00080 


c 


■ ■ WILL BE PLOTTED ON GRID NO,l 

GRFODO^O 


c 


TF L IS POSITIVE - NP POINTS WILL BE PLOTTED ON GP|0 2 riR 

BGftfOOlOD 


c 


IF L IS ZERO - NP POINTS WILL BE PLOTTED ON PREVIOUS GRID 

GftFOOllO 


c 



GPF0012D 


c 


ISYM PLOTTING SYMBOL 

GRF00130 


c 



GRF00140 




fiCDX ALPHA-NUMERIC CHARACTERS FOR THE X AXIS L Y AXIS.. THESE 

GRFD015D 


c 


BCOY ARRAYS ARE LIMITED TO 72 CHARACTERS EACH AND THE LAST 

GRFODtAO 


c 



GRF00170 


c 


NP NUMBER OF POINTS TO BE PLOTTED 

GRF00180 


c 


IF NP IS negative - THE POINTS WILL BE CONNECTED, THUS THE 

GPFDD19D 


c 


VALUES OF X MUST BE IN ASCENDING ORDER 

GRF00200 


c 



GRF002L0 


c 


XMAX MAXIMUM C MINIMUM VALUES OF X C Y, TF BOTH MINIMUM AND 

GRF00220 


c 


XMIN MAXIMUM VALUES OF X ARE ZERO, THE ARRAY IS SEARCHED FOR 

GRF00230 


c 


VMAX THE APPROPRIATE VALUES* THE SAME PROCEDURE IS FOLLOWED 

GftP002A0 


c 


YMIN FOR MINIMUM £ MAXIMUM Y 

GRF00250 


c 



GftF00260 

000? 



DIMENSION XI5001 , Y( 50DJ ,BCOX( 181 tBCDYdSl 

GPF0D2‘'0 

0003 . 



NZ * lABSINP) 

GRF00275 

000^ 



IF(NZ .FO* 01 GO TO 60 

GRF00280 

000'5 



IFU *E0* 0 1 GO TO 30 

GRF0029P 

0006 



XL = XMIN 

GRF00300 

0007 



Xft * XMAX 

GRF00310 

0005 



IFUXL *NE. O.OI .OK. (XR ,NE* 0*011 GO TO 15 

GRF00320 

000^ 



XL = X( 1) 

GPF00325 

0010 



XR = xm 

GRF00326 

0011 



DO 10 K i l,NZ 

GRF0033D 

0012 



XL * AMINH X( K1 ,XL1 

GRF00340 

0013 


10 

XR = AMAXK XIKJ.XRI 

GRFO035D 

0014 


15 

Yfl = YMIN 

GP^00360 

0015 



YT = YMAX 

GftF003’^0 

001 6 



IFMVT .NE* O.OI .Oft. (YB .NE. O.OM GO TO 25 

GftF 00380 

0017 



YB = VU I 

Gftf00385 

0018 



YT = Y(l) 

GRF003R6 

0019 



on 20 J » ItNZ 

GRF004CO 

0020 



yp « AMI N1 1 YUl tYBI 

GPF00410 

0021 


20 

YT = AMAXl ( V( J1 tYTI 

GftF00420 

0022 


25 

IFIL .LT* 0 1 LI * 1 

GftF00440 

0023 



TFfL .GT. 0 1 LI = 2 

GRF00450 

0024 



NCX = N8LANC(0CDXtlfl) 

GftF00460 

0025 



NCY = NBLANCOCDV.ISI 

GRF 00470 

0026 



LGO » lABSlLP 

GRFOO40O 

0027 



OCX » 10*0 

GRF00490 

0028 



ncv * 10. 0 

GRF 00500 

0029 



INCRY = -14 

GRF 00510 

0030 



call MARGNR IL* ICY) 

GRF00520 

0031 



IX « 524 - 44NCX 

GPF00530 

0032 



lY = ICY 4T4NCY 

GPF00540 

0033 



IFILGO ,E0. 1 1 lY.l = 0 

GRF00550 

0034 



TF(LGO *E0* 21 lYl * ICY - 253 

CftF00560 

0035 



TFUGO -FO* 31 lYl * ICY - 169 

GPF00570 

0036 



CALL OXDYVf 1 , XL , XR t DX tN, 1 , NX * DC X , I ER I 

GRF 00 580 

0037 



CALL OXDYV(2,YBfYTfOYtMt J,NY,OCY, lERI 

GRF00590 

0038 



CALL GRTniV(Ll,XLtKft,yB,YTtOX*OY*NtM,I, JtNX.NYl 

GRF00600 

0039 



CALL PRTNTVlNCXtBCnXtlXtlYll 

GRF0n6l0 

0040 



CALL APRNTVtO, INCRV,NCVt0COY,OtlYl 

GRF00620 

0041 


, 30 

CALL APLflTVfNZ,XtYtltltlttSYMtIER> 

GRF 00640 

0042 



IF <NP .GT, 01 GO TO 60 

GRF 00650 

0043 



DO 50 KK = IfNZ 

GRF 00^60 

0044 



NXA = NXV(XIKK) 1 

GRF00670 

0045 



NVA = NYVIYCKKM 

GRF006AO 

0046 



IF(KK *LT, 2) GO TO 40 

GRF006«50 

0047 



CALL LlNEVINXOtNYOfNXA.NYAI 

GRF 00 700 

0048 


40 

NXO » NXA 

GRFOOYIO 

0049 



NVO - NYA 

GRF00'^?0 

0050 


50 

CONTINUE 

GRF0D730 

0051 


60 

RETURN 

GRF00740 

0052 



END 

GRF00750 
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000 1 


0002 

0003 

0004 

0005 

0006 
C007 
0008 

0009 

0010 
oou 
0012 


c 

c 

c 


Sll8R0UTIMe MARGMRfLltMtJI MARGOOlO 

HAftG 0020 

‘^ MAftGNR KARGINS^FGR ■pLOn 3 GRIDS PEP FRAME MARG0030 

MARG^040 


DlMEMStON ■U6,3>, HARG0050 

DATA L /24f 35Tf 690,24 t524f 24, 600,357,34,524^24* 24,1 T8 ,511 t844t262tHAPG0060 


I 762,512/, IN /5f3,0/ 

_ 1F(__LI *GF, 0 I GO TO 10 

-Cl 

N » 1 

10 rc = N + IN(K) 

" ^ ML3 = Li iC,31 

call SETMIVI24,0,UTCTU,L(tC,2n 
N 

R^tUQN“ ' ' ' 

END 


MART, 0070 
MAPG0080 
MARGOOOO 
MARGOlOO 
MARGOllO 
-MARG0120 
MARG0130 
MARGO 140 
MARGO 150 
MARC0160 


0001 



FUNCTION NBLANC(WORD,NI 

NBLKOOlO 


C 



NBLK0020 




•NRLANC* DETERMINES NO. OF ’CHARACTER S IN A MOLLERITH LABEL 

NBLK0030 


C 


IF THERE ARE NO CHARACTERS, NBlANC ^ 0 

NRLK0040 


C 



NBLK0060 

0002 



DIMENSION NOROtlftl 

NBLKOOTO 

0003 



DATA BLANK /• '/ 

N8LKOORO 

0004 



no 10 M =< UN 

NBLK0090 

ooo's 



i^ = N - H 

NBLKOlOn 

0006 



iFMMORDf I >-dLANKl .NE* 0.0» GO TO 20 

N6LK0110 

0007 


10 

CONTINUE 

NBLK0X20 

0008 



I = 0 

N0LKO13O 

0009 


20 

N6LANC = 441 

NRLK0140 

0010 



RETURN 

NBLK0150 

OOU 



'end' ■ '■ ■ ^ 

N0LKO16O 
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